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Abstract. We start with a review of a class of systems with invariant relations, so 
called systems of Hess-Appel'rot type that generalizes the classical Hess-Appel'rot 
rigid body case. The systems of Hess— Appel'rot type have remarkable property: there 
exists a pair of compatible Poisson structures, such that a system is certain Hamil- 
tonian perturbation of an integrable bi-Hamiltonian system. The invariant relations 
are Casimir functions of the second structure. The systems of Hess-Appel'rot type 
carry an interesting combination of both integrable and non-integrable properties. 

Further, following integrable line, we study partial reductions and systems having 
what we call the Zhukovskii property: These are Hamiltonian systems on a symplectic 
manifold M with actions of two groups G and K; the systems are assumed to be K— 
invariant and to have invariant relation <E> = given by the momentum mapping of 
the G-action, admitting two type of reductions, a reduction to the Poisson manifold 
P = M/K and a partial reduction to the symplectic manifold No = < I >— 1 (0)/G?; 
final and crucial assumption is that the partially reduced system to No is completely 
integrable. We prove that the Zhukovskii property is a quite general characteristic of 
systems of Hess-Appel'rote type. The partial reduction neglects the most interesting 
and challenging part of the dynamics of the systems of Hess-Appel'rot type - the 
non-integrable part, some analysis of which may be seen as a reconstruction problem. 

We show that an integrable system, the magnetic pendulum on the oriented 
Grassmannian Gr + (n, 2) has a natural interpretation within Zhukovskii property and 
that it is equivalent to a partial reduction of certain system of Hess-Appel'rot type. 
We perform a classical and algebro-geometric integration of the system in dimension 
four, as an example of a known isoholomorphic system - the Lagrange bitop. 

The paper presents a lot of examples of systems of Hess-Appel'rot type, giving 
an additional argument in favor of further study of this class of systems. 
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1. Introduction 

Historically, the Hess-Appel'rot system as a classical rigid-body system, appeared 
(see [35]) just a year after the celebrated Kowalevski 1889 paper [38], and its imme- 
diate popularity had been connected with its relationship with the Kowalevski paper. 
Kowalevski started, as we know, from a careful analysis of the solutions of the Euler 
and the Lagrange case of rigid-body motion and formulated a problem of describing 
the parameters (A, B, G, xq, y^, zq), for which the Euler-Poisson equations have a general 
solution in a form of uniform functions only with moving poles as singularities. Here, 
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/ = diag(A, B, C) represents the inertia operator, and \ = { x OiVOi z o) is the centre of 
mass of the rigid body. 

Then, in §1 of [38], after some necessary conditions were formulated, she descovered 
a new case, which is now known as the Kowalevski case. The last case was, according 
to Kowalevski, a unique possible beside the cases of Euler and Lagrange. However, 
considering the situation where all momenta of inertia are different, Kowalevski came to 
a relation analogue to the following (see [34]): 

xoVMB -C) + yoVB(C -A) + z^C(A -B) = 0, 

and she concluded that xq = yo = zq, which represented the Euler case. 

However, Appel'rot noticed few years later, that the last relation admitted one more 
case, not observed in [38]: 

x oy /A(B -C)+ z ^/C(A -B) = 0, 2/o = 0, 

where he assumed A > B > C. Such intriguing position corresponding to the possible 
mistake in the Kowalevski paper, made the Hess- Appel'rot systems very attractive for 
leading Russian mathematicians from the end of XIX century as a possible counterex- 
ample. But, after a few years, Nekrasov and Lyapunov proved that the Hess- Appel'rot 
systems didn't satisfy the condition investigated by Kowalevski, which means that con- 
clusion of §1 of [38] was correct. 

And, from that moment until very recently, the Hess- Appel'rot systems were basically 
left aside, even in modern times, when new methods of inverse problems, Lax represen- 
tations, finite-zone integrations were applied to almost all known classical systems. 

A modern theory of systems of Hess- Appel'rot type has been developed in [18, 19, 
20, 21]. It started with a construction of a Lax representation for the Hess- Appel'rot 
system in [18], see the Proposition 2.1 below. Generalization of this Lax pair in four- 
dimensional case led to construction of a new integrable rigid-body system in [18], called 
the Lagrange bitop. Algebro-geometric integration procedure of the Lagrange bitop has 
been performed in [19]. It brought to the discovery of a new class of integrable systems, 
which was named isoholonomic systems in [19]. Higher dimensional generalizations of 
the classical rigid-body Hess-Appel'rot systems have been constructed in [20, 21] as cer- 
tain perturbations of the isoholonomic integrable systems. Finally, in [21] after detailed 
analysis of infinite set of new examples in arbitrary dimensions, the theory of systems of 
Hess-Appel'rot type has been settled down providing an axiomatic, general and abstract 
approach, see also the Section 3. 

According to this theory the systems of Hess-Appel'rot type form a class of dynamical 
systems, obtained as certain perturbations of integrable, bi-Hamiltonian systems which 
carry an interesting combination of both integrable and non-integrable properties. 

Suppose a bi-Poisson structure {•, -}i + A{-, -}2 is given, with an integrable, bihamil- 
tonian system with the Hamiltonian H corresponding to the first structure. Further, let 
/i, . . . , fk be the commuting integrals of the system (Ho, {•, - }i), which are Casimirs for 
the second structure {•, -}2- Then, a system of Hess-Appel'rot type is Hamiltonian with 
respect to the first structure with a Hamiltonian 

k 

H = H + Y,Jibifi, 
i=i 

where J; are constants and 6; are certain functions on the phase space. The invariant 
relations are 

fi = 0, l = l,...,k. 
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Thus, the invariant manifolds are symplectic leaves of the second Poisson structure. 

As perturbations, they are global and not just small perturbations as it is usually 
the case in the study of non-integrable perturbations of intcgrable systems. A balance 
between integrable and nonintegrable properties is obtained by the choice of perturba- 
tions. The system of Hess-Appel'rot type is Hamiltonian with respect to the first Poisson 
structure, but perturbations and invariant relations are defined by Casimirs of the second 
structure. 

For the classical Hess-Appel'rot case, its integrable part made it close to the Kowale- 
vski study while its nonintegrable side finally disqualified it as a possible counterexample 
for the Kowalevski statement. Classical integration of its integrable part one may find also 
in the book of Golubcv [34] , while algcbro-gcometric integration has been performed in 
[18]. ft was Zhukovskii who observed (see [73]), that after certain reduction, the classical 
Hess-Appel'rot system reduces to the completely intcgrable system of spherical pendulum 
(see the Subsection 5.1 below). This observation of Zhukovskii motivated us to introduce 
the notion of Zhukovskii property, see the Section see Section 5. 

For the four-dimensional generalization of the Hess-Appel'rot system, detailed sep- 
aration of integrable and nonintegrable part has been done in two ways in [20, 21], 
both classically and algebro-geometrically. Moreover, two integration procedures for the 
integrable part have been performed in all details in [21]. 

From these integrations, one can see that a completely integrable system in the same 
space requires one integration more in three-dimensional case and in the four-dimensional 
case it would require two integrations more. 

Reviews of these results can be found in [17] and [30]. 

Following [21], the interest for the Zhukovskii property in the modern context has 
been expressed in [36, 37], where a study of partial reductions has been developed further, 
see Section 5. 

One of the aims of the present paper is to provide a systematic study of the Zhukovski 
property from a general point of view of the axioms of the systems of Hess-Appel'rot 
type, ft appears that the Zhukovskii property is a quite general characteristic of the 
systems of Hess-Appel'rot type and together with the partial reduction it traces well the 
integrable part of a system of Hess-Appel'rot type, see Sections 5, 6. In the same time, it 
neglects totally the most interesting and challenging part of the dynamics of the systems 
of Hess-Appel'rot type - the non-integrable part. This blindness to the non-integrable 
part, makes the partial reduction being important but of a limited range and domain in 
complete understanding and studying of the systems of Hess-Appel'rot type. 

Nevertheless, a completely integrable system, the magnetic pendulum on GV + (4, 2) 
which has been introduced and studied in [14] within the study of magnetic flows on 
homogeneous spaces, appears to have an interpretation within Zhukovski property: we 
show that it can be obtained as a partial reduction of a certain system of Hess-Appel'rot 
type. Moreover, it appears that the magnetic pendulum is equivalent to a very simple 
instant of the Lagrange bitop from [19]. Thus, the magnetic pendulum is an example of 
an isoholonomic system. So, the integration techniques of [19] and [21] can bee applied 
directly to the magnetic pendulum. Complete integration of the magnetic pendulum in 
classical and algebro-geometric manner is the second main aim of this paper and it is 
performed in the Section 8. 

The paper provides a lot of examples of systems of Hess-Appel'rot type, motivated 
from [20, 21, 36, 37, 23] and from references therein, see the Sections 4 and 7. Such 
a rich set of examples is additional argument in favor of further study of the systems of 
Hess-Appel'rot type. 
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At the end of the paper, we collected for a reader's sake all necessary notions from 
the theory of Hamiltonian systems and their reductions, see the Appendix. 

2. Classical three-dimensional Hess-Appel'rot system 

The Euler-Poisson equations of the motion of a heavy rigid body in the moving frame 
are [34]: 

M = MxJ7 + Txx, 

(2.1) t = Txfl 

n = JM, J = diag(Ji, J 2 , J 3 ), 

where M is the kinetic momentum vector, $7 the angular velocity, J a diagonal matrix, 
the inverse of inertia operator, T a unit vector fixed in the space and \ is the radius 
vector of the centre of masses. 

It is well known (see for example [34]) that equations (2.1) have three integrals of 
motion: 

(2.2) F 1 = i(M,fi) + (r, x ), F 2 = (M,r), F 3 = (r,r) = i. 

Thus, for complete integrability, one integral more is necessary [34] . Let Ji < J 2 < J3 
and x — (^07 2Ah ^o). Hess in [35] and Appel'rot in [4] found that if the inertia momenta 
and the radius vector of the centre of masses satisfy the conditions 

2/o = 

(2-3) rz rz r 

xq v J3 - h + z y J 2 - Ji = 0, 

then the surface 

(2.4) Fi = M lXa + M 3 z = 

is invariant. 

The compact connected components of the regular invariant sets given by (2.2), 
(2.4) are tori, but not with quasi-periodic dynamics. The classical and algebro-geometric 
integration can be found in [34] and [18], respectively. It is shown that the equations of 
the motion reduce to one elliptic integral and one Riccati differential equation. 

The Zhukovskii geometric interpretation of the conditions (2.3) [73, 40]. Let us 
consider the ellipsoid 

M? Ml Ml 
-J- + -J- + -J- = 1, 

J\ J2 J3 

and the plane containing the middle axis and intersecting the ellipsoid at a circle. Denote 
by / the normal to the plane, which passes through the fixed point O. Then the condition 

(2.3) means that the centre of masses lies on the line 

Having this interpretation in mind, we choose a basis of moving frame such that the 
third axis is /, the second one is directed along the middle axis of the ellipsoid, and the 
first one is chosen according to the orientation of the orthogonal frame. In this basis (sec 
[15]), the particular integral (2.4) becomes 



F 4 = M 3 = 0, 
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and the matrix J and mass centre \ obtain the form: 

/ Jl J13 \ 

(2.5) J= Jr , x=(0,0,z ). 

V ^3 o j 3 y 

This will serve us as a motivation for a definition of the four-dimensional Hess-Appel'rot 
system. 

A three-dimensional Lagrange top is defined by the Hamiltonian: 

1 / M? + Mo M|\ 

according to the standard Poisson structure 

{Mi,M,}i = -e ijk M k , {Mi, Tj}i = -e ijk T k , {T^Tj} = 

on the Lie algebra e(3). It is also well-known that three-dimensional Lagrange top is 
Hamiltonian in another Poisson structure, compatible with first one. This structure is 
defined by: 

{r^r^h = -e ijk T k , {M ll M 2 } 2 = 1, 
and the corresponding Hamiltonian is: 

H L = (a- 1)M 3 Q(M 2 + Ml) + L 3 ) + M 1 T 1 + M 2 T 2 + M 3 T 3 

where ii = 1, J 3 = a, z = 1. Casimir functions in the second structure are Lf + T 2 + P 3 
and M 3 . 

Let us observe that the Hamiltonian for the three-dimensional Hess-Appel'rot case 
is a quadratic deformation of Hamiltonian Hl of the Lagrange top: 

Hha = H l + J 13 MiM 3 . 

The function M 3 , which gives the invariant relation for the Hess-Appel'rot case, is a 
Casimir function of the second Poisson structure. 

Let us mention that the Lax representation for Hess-Appel'rot system (2.1), (2.3), 
(2.4) is constructed in [18]: 

Proposition 2.1. ([18]) On invariant manifold given by the invariant relation (2.4), the 
equations of Hess-Appel 'rot system are equivalent to the matrix equation 

L(X) = [L(X),A(X)} 

where L(X) = X 2 C + AM + T, A(X) = X X + n. 

Here we denoted with M £ so(3) antisymmetric matrix that corresponds to the 
vector Mel 3 due to correspondence Mjj = —tij k M k (and similar for C, T, 0,%), and 
C = -j~X- Another Lax representation for the Hess-Appel'rot system is given in [22]. 
Using it a sort of separation of variables for Hess-Appel'rot system is discussed there. By 
putting n = 3 in (4.12) we get yet another Lax representation (see [36]). 

3. General systems of Hess-Appel'rot type. Axiomatic approach 

3.1. The first set of axioms: general Poisson settings. Suppose a Poisson 
manifold (M 2 ™, {•,•}) is given, together with k + 1 functions H,f u ...,f k e C°°(M), 
such that 
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(Al) 

k 

= 8,6^(1), i,j = l,...,k; 

i=i 

(A2) 

{/»/?} = 0, i,j = l,...,k. 
A more general case can be obtained by replacing condition (A2) with 

(A2') 

k 

fj} = 4*^' 4* = C ° nSt > = lj k - 

2 = 1 

In this case, the algebra of invariant relations is a noncommutative Lie 
algebra. 

Starting from the Hamiltonian system (M, Hq) with k integrals in involution 
fi,...,fk, choosing functions bj G C°°(M), j = l,...,k, one comes to a 
rcstrictively integrable system: 
(HP) (Hamiltonian perturbation axiom) 

The system (M, H) where 

k 

H = H + J2b j f j , 
i=i 

will be called a Hamiltonian perturbation. It satisfies (Al) with 

a ij = {Pji fi}i *i J = 1) • • • > k. 

(BP) (Bi-Poisson axiom) There exist a pair of compatible Poisson structures, such 
that the system is Hamiltonian with respect to the first structure, having the 
Hamiltonian of the form (HP), such that fi are Casimir functions with respect 
to the second structure. 

The invariant relations define symplectic leaves with respect to the second structure, 
and the system is Hamiltonian with respect the first one. 

3.2. The second set of axioms: Kowalevski property. To get the right choice 
of axioms, we have to turn back to the Kowalevski analysis. First, we are going to 
introduce some general notions, see [39]. 

Suppose a system of ODEs of the form 

(3.1) Zi = fi(zi, . . . ,z n ), i = l,...,n, 

is given and there exist positive integers gi, i = 1, . . . , n, such that 

/i(a Sl zi, . . . ,a 9n z n ) = a 9 ' +1 fi(zi, z n ), i = l,...,n. 

Then the system (3.1) is quasi-homogeneous and numbers gi are exponents of quasi- 
homogeneity. Then, for any complex solution C = (C\, . . . , C n ) of the system of algebraic 
equations: 

(3.2) -g i C i = fi(C 1 ,...,C n ), i = l,...,n, 
one can define the Kowalevski matrix K = K(C) = [K % -{C)\. 
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Eigcn- values of the Kowalevski matrix are called the Kowalevski exponents. This ter- 
minology was introduced in [71]. In last twenty years, heuristic and theoretical methods 
in application of Kowalevski matrix and Kowalevski exponents in study of integrability 
and nonintegrability have been actively developing, see for example [2, 3, 39]. But, the 
notion of Kowalevski matrix and Kowalevski exponents were introduced by Kowalevski 
herself in [38]. The criterion she used (see [38], p. 183, 1. 15-22) to detect a system 
which is now known as the Kowalevski top, can be formulated in Yoshida terminology 
as: 

Kowalevski condition (Kc) . The 6x6 Kowalevski matrix should have five different 
positive integer Kowalevski exponents. 

Now we return to the study of systems of Hess-Appel'rot type. The systems we have 
constructed are quasi-homogeneous. Exponents of each M variable are g = 1, and for 
any T they are equal to two. We are going now to calculate Kowalevski exponents for 
the Hess-Appel'rot systems. 

Three-dimensional Hess-Appel'rot case. Let us denote (-Mi, M3, T\, I^, 
T3) by (zi, . . . , z@). Then the Euler-Poisson equations take the form (3.1) with 

fi = (J3 - Ji)z 2 z 3 + J13Z1Z2 + z 5 ; 

h = (J3 - Ji)ziz 3 + Ji 3 (z% - z\) - z 4 ; 

h = J13Z2Z3; 

fi = J3Z5Z3 - JlZ 2 Z e + J13Z1Z5; 

h = -J3Z3Z4 + J\Z\z e + Ji 3 (z 3 z 6 - Z1Z4); 
fe = J3Z2Z4. — J\Z\z b - J13Z3Z5; 

and <7j = 1, i = 1,2,3 and g$ = 2, i = 4,5,6. The invariant relation corresponds 
to the constraint C3 — 0. So, we are looking for solutions (ci, C2, 0, C4, C5, c%) of the 
system of the form (3.2). One can easily get C4 = —J\3c\ + c 2 , c 5 = — ci(l + J13C2), 
cq = — (cf + c\)/2. Then, for c\ 7^ 0, we get four possible solutions for (ci,C2) divided 
into two pairs: (±z/Ji 3 , — I/J13) and (±2?/Ji3, — 2/J13). The Kowalevski exponents are 

(-1,-2,2,4,3,3), (-1,1,3,2,2,2), 

respectively. 

Thus, it can easily be seen that classical Hess-Appel'rot system doesn't satisfy exactly 
the Kowalevski condition (Kc), although it is quite close to. 

Thus, using into account properties of Kowalevski exponents of algebraically- inte- 
grablc Hamiltonian systems, we can conclude that for the systems we have constructed, 
functions hi in the perturbation formula (HP) should satisfy two conditions: 

(QH) (quasi-homogeneity axiom) The obtained system of Hamiltonian equations 
has to be quasi-homogeneous. 

In such Kowalevski matrix exists and wc come to the last condition. 

Suppose the invariant relations correspond to equations z\ = 0, . . . , z& = 0. 

Denote by p number of Casimirs: n = p + 2m, where 2m is the dimension 
of a general symplectic leaf. 
(ArA) (Arithmetic axiom) For any nonzero solution C — (0, 0, Ck+i, c n ) of the 
system (3.2), the Kowalevski matrix K(C) has n—p eigen-vectors tangent to the 
symplectic leaf and p transversal to it. Half of the Kowalevski exponents which 
correspond to tangential eigen-vectors and all of transversal ones are rational 
numbers. Irrational numbers among the second half of tangential Kowalevski 
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exponents are divided into pairs such that the differences are integrally depen- 
dent. 

The axioms of systems of Hess-Appel'rot type provide conditions which determine 
classical Hcss-Appel'rot system among three-dimensional systems of Hess-Appel'rot type. 
More precisely, suppose the two Poisson brackets are given on e(3) as above and a system 
is given by a Hamiltonian 

(3.3) H 1 =H + JbM 3l 

where Ho is the Hamiltonian of the Lagrange top corresponding to the first Poisson 
structure, M 3 is its integral and a Casimir for the second structure, J is a nonzero 
constant and & is a function, such that the axioms of the systems of Hess-Appel'rot are 
satisfied. 

In [21] the following rigidity theorem has been proved. 

Theorem 3.1 ([21]). The only non-zero polynomials b which give systems of Hess- 
Appel'rot type by relation (3.3) are of the form 

b(zi, . . . ,z 6 ) = zi + kz 3 . 

All systems of Hess-Appel'rot type of the form (3.3) are the classical Hess-Appel'rot sys- 
tems. 

The last theorem provides a strong argument in favor of the choice of axioms which 
have been postulated in [21]. We take these axioms as the starting point of our current 
research. 

Since the Zhukovskii property is main object of our study in this paper, we are going 
to focus ourselves on the first set of axioms. 

4. Examples of systems of Hess-Appel'rot type 

4.1. Rigid body systems on so(n) x so(n). The Euler-Poisson equations of motion 
of a heavy rigid body fixed at a point are Hamiltonian on the Lie algebra e(3), which is 
the semi-direct product of Lie algebras R 3 and so(3). Since M 3 is isomorphic to so(3), 
there are two natural higher-dimensional generalizations of Euler-Poisson equations. The 
first one is given by Ratiu in [64] and it is to the semi-direct product so(n) x so(n) and 
the second one is to the Lie algebra e(n) = M™ x so{n). 

Equations of a heavy n-dimensional rigid body on so{n) x so(n) are : 

M=[M,n] + [T,x] 

(4.1) 

r = [r,o], 

where M, f2,L, x £ so(n), il = AM and x is a constant matrix (sec [64]). Here A : 
so(n) — > so(n) is the inverse of the rigid body kinetic energy operator. The Euler-Poisson 
equations (4.1) are Hamiltonian with the Hamiltonian function 

(4.2) H = l(M,fi> + ( X ,r) = -±tr(Mfi) - itr( X L), 

in the standard Poisson structure on the semi-direct product so(n) x so(n): 

{M ij ,M jk } 1 = -M ik , {M lJ ,r Jk } 1 = -v ik , {r y ,r fcJ }i = 0. 

The Casimir functions are tr(r 2fc ), tr(MT 2fc+1 ) and the dimension of generic sym- 
plectic leaf is n(n — 1) — 2 [^] . 
We will suppose that 

(4.3) n = JM + MJ, 
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where J is a constant symmetric matrix. The operator M i— ► JM + MJ belongs to the 
class of Manakov operators [50] on so(n). 

The Lax representation together with Zhukovskii geometric interpretation presented 
in the second section, where inspiration for a construction of a higher-dimensional gen- 
eralization of Hess-Appel'rot system in [21]. Let us first consider the four-dimensional 
case. 

Definition 4.1. The four-dimensional Hess-Appel'rot system is described by the equa- 
tions (4.1), (4.3) and conditions: 



(4-4) 



J 



( Ji 
Ji 

Jl3 

V o j 24 

such that the operator M 





J 3 





J24 


■h) 



( 

-X12 



V 



X12 







\ 


X34 
-X34 j 



JM + MJ is positive definite and X12 + X34 0- 

The invariant surfaces are determined in the next lemma. 

Lemma 4.1. (i) For the four- dimensional Hess-Appel'rot system, the following relations 
take place: 

M12 = Ji 3 (M 13 M 12 + M 2i M 3i ) + J 2i (M 13 M 3i + M 12 M 24 ), 

M 34 = Ji 3 (-M 13 M 34 - M l2 M 2A ) + J 2 ±(-M 13 M 12 - M 2A M 3A ). 

(ii) The system has two invariant relations: 

(4.5) M12 = 0, M 34 = 0. 

Now let us introduce a new Poisson structure, compatible with the standard one, as 
follows: 

{r«,r J - fc } 2 = -r ifc , {M y ,r fcJ } 2 = o, {m 13 ,m 23 } 2 = - X i 2 , 

{M14, M 24 } 2 = -X12, {M13, Mi 4 } 2 = -X34, {M 23 , M 24 } 2 = -X34- 
Casimir functions in this structure are M12, M 34 , T\ 2 + Tf 3 + Y\ A + T 23 + T| 4 + r 34 , and 

Ti 2 r 34 -i- r 23 ri 4 — ri 3 r 24 . 

The situation with four-dimensional Hess-Appel'rot case is similar to the three-dime- 
nsional case: the Hamiltonian is again a quadratic deformation: 

Hha = H LB + Ji 3 (-M 12 M 23 + Mi 4 M 34 ) + J 24 (Mi 2 Mi 4 - M 23 M 34 ), 

where H LB = \ (2 J X M? 2 + (Ji + J 3 )M 1 2 3 + (Ji + J 3 )M? 4 + (Ji + J 3 )M| 3 + (Ji + J 3 )M| 4 + 
2J 3 M| 4 ) + Xi2^i2 + X34T34 is the Hamiltonian function of the Lagrange bitop. The 
Lagrange bitop is a complete intcgrable system of a heavy rigid body on so(4) x so(4) 
defined in [18] and studied in details in [19]. Moreover, the Lagrange bitop is a bi- 
Hamiltonian system. Assume that J\ — a, J 3 = 1 — a. Then 

(2a- l)(xi 2 M 12 +x 34 M 34 ) 



(4.6) 



Hlb 



X12 



X34 



M| 3 



M| 4 



+ Xl2Tl2 + X34r 34 



+ 



(l-2a)(xi2M 34 + X34Mi 2 ) 



(M 23 Mi 4 - Mi 3 M 24 + xi2r 34 + X34r i2 ) 

+ M 12 r 12 + M 13 r 13 + m 14 l 14 + M 23 r 23 + m 24 l 24 + M 34 r 34 . 



x?2 



X34 



is the Hamiltonian in the second structure (4.6). The functions M 12 and M 34 , giving 
invariant relations for the four-dimensional Hess-Appel'rot system, are Casimir functions 
for the Poisson structure (4.6). 



SYSTEMS OF HESS-APPEL'ROT TYPE AND ZHUKOVSKII PROPERTY 



11 



Let us pass to the arbitrary dimension n > 4. 

Definition 4.2. The n-dimensional Hess-Appel'rot system is described by the equations 

(4.1), (4.3), together with conditions: 

(4.7) 

J = diag(Ji, Ji, J 3 , . . . , J 3 ) + Ji 3 (Ei <Z>E 3 + E 3 <Z> E{) + J 24 (E 2 <Z>E 4 + E 4 <Z> E 2 ), 
X = X12E1 AE 2 , Xi2 7^0. 

Invariant relations are given in next lemma. 

Lemma 4.2. (i) For the n-dimensional Hess-Appel'rot system we have: 

n 

Mil =Ji 3 (M 12 M 13 + M 24 M 34 + M 2p M 3p )+ 

p— 5 
n 

J 24 (M 12 M 24 + M 13 M 34 - M ip M 4 P ) 

p— 5 

71 

M34 = - Ji 3 {M 13 M 34 + M 24 M 12 + £ M lp M p4 )- 

p—5 

n 

J 24 (M 13 M 12 + M 24 M 34 + M 2p M 3p ) 

p— 5 

M 3p = - Ji 3 {M 13 M 3p + M 2p M 12 ) - J 24 (M 34 M 2p + M 23 M 4p )+ 



M 34 n 4p - n 34 M 4p + ^2(M 3k n kp - n 3k M 4p ), P > 4, 

fc=5 

M 4p =Ji 3 (-M 14 M 3p + M lp M 34 ) + J 24 (M 12 M lp - M 24 M 4p )- 

n 

M 34 n 3p + n 34 M 3p + J2(M 4k n kp - n 4k M 4p ), P > 4, 



k=5 

Mm =0, k, I > 4 



(ii) The system has the following set of invariant relations: 

(4.8) M 12 = 0, M lp = 0, l,p>3. 

In [64] the n-dimensional Lagrange top on the semidirect product so(n) x so(n) is 
constructed. In the metric = JM + MJ, where J = diag(Ji, Ji, J 3 , . . . , J 3 ), the 
n-dimensional Lagrange top is defined with Hamiltonian 

Hl = \[ 2JiM* 2 + (.h + j 3 )J2( m i p + M D + 2 ^3 £ Ml q + y 12 r 12 

y p=3 3^p<q^n J 

Consider a Poisson structure 

(4.9) {T ij ,r jk } 2 = -r ik ,{M ij ,M kl } 2 = 0,{M ll ,M 2l } 2 = -l, i = 3,...,n, 
compatible with the standard one. The dimension of a symplcctic leaf in this structure 



is 



( "- 2)( "~ 3) - [2=2] + 4(n - 2), hence there arc «!=|™+8 + [«] Casimir functions: 



Mi 2 , M M , tr(r 2fc ), 2<p<q<n, k = 1, . . . , 
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The n-dimensional Lagrange top is also bi-Hamiltonian system. In the Poisson struc- 
ture (4.9) its Hamiltonian is: 




H L = (2a - 1)M 12 - ±JMi p + Mi p ) + T 12 



[l -2a) ^2 M pq (M lq M 2p - M 2q M lp + T pq ) + ^2 M pq T pq, 
3^p<q^n l^.p<q^.n 

where J\ = a, J 3 = 1 — a, \12 = 1- 

Similarly as in dimension 3 and 4, Hamiltonian for the Hess-Appel'rot system in 
arbitrary dimension n is a quadratic deformation of the Hamiltonian for the n-dimensional 
Lagrange top: 

n 

Hra = H l + ^2(J 13 M lk M 3k + J 24 M 2k M 4k ), 
k=l 

and functions M\ 2 , M pq ,p, q ^ 3, which give the invariant relations, are Casimir functions 
for the Poisson structure (4.9) 

Next theorem gives a Lax pair for the Hess-Appel'rot system. 

Theorem 4.1 ([21]). On invariant manifold given by the invariant relations, the equa- 
tions of n-dimensional Hess-Appel'rot system are equivalent to the matrix equation 

L(X) = [L(X),A(X)}, L(\)=\ 2 —^— x + \M + T : A{\)=\ x + n. 

J\ + J3 

Using this Lax representation, the both classical and algebro-geometric integration 
procedures are presented in [21] in dimension four. 

4.2. Rigid body systems on e(n). Let us now consider rigid body motion on the 
Lie algebra e(n). The standard Poisson structure on e(n)* = e(n)(M,T) is given with: 

{Mij, M w }i = -5 jk Mu, Tj}! = 0, {My ,r fc }i = -Ti6 jk + Yj5 ik . 

Here we identified e(n)* = e(n) by the use of a non-invariant scalar product: 

1 " 

<(M, r), (m, r)> = (m, m) + (r, r) = --tv(MM) + 

1 i=l 

A heavy rigid body Hamiltonian read: H = i(M, O) + (x,T), where \ G ^™ is the 
vector of mass centre and T G R n is a vertical vector considered in the moving coordinate 
system and, as above, = AM. We can choose x = XnE n — (0, . . . , 0, \ n ). The 
corresponding Euler-Poisson equations are: 

M = [M, Q] + XnE n A r 

(4.10) 

r = -or. 

In [7] Belyaev considered the n-dimensional Lagrange top defined by the Hamiltonian 
function: 

(4.11) H A = ^(M a ,M D ) + ^(M Sl M 3 )+ Xn r n , 

where so(n) — g © is orthogonal symmetric-pair decomposition of so(n): g = (Ei A 
Ej 1 1 < i < j < n — 1) = so(n — 1), = (Ei A E n 1 1 < i < n — 1) and a, b > are real 
parameters. Note that the kinetic energy has the Manakov form (4.3), where we take 
= J A M + MJ A , J a = diag(Ji, J 1 ,...,J 1 , J n ). Then a = J x + J n , b = 2J X . 
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Bclyaev proved noncommutative integrability of the system [7]. The Lax represen- 
tation is given by Reyman and Semenov-Tian-Shanski [66] : 

(4.12) L(X) = [L(X),A(X)], L(X) = f + AM + X 2 ^-E n , A(X) = lo + X X E n . 
Here, for a given (M, T) e e(n), M,T G so(n + 1) are defined by 

Let £ be the set of integrals obtained from (4.12) 

(4.13) £: tr(L(X) 2k ), k = 1, . . . , rank SO(n + 1), AeK 
and let S be linear functions ong = so(n — 1) 

(4.14) S = {M i3 1 1 < i < j < n - 1}. 

Then {£,£}i = 0, {£, S}\ = and £ + S is complete set of integrals of the Lagrange 
top system (4.10), (4.11). 

As above, let us define a Hess-Appel'rot system on e(n) as a perturbation of the 
Lagrange top using operator (4.3) with 

(4.15) J = J A + Jn, Ja = diag(Ji, Ji,..., Ji, J„), Jn = Jin{E\ ® E n + E n ® E\). 
Proposition 4.1. TTie equations (4.10), (4.3), (4.15) have invariant relations: 

(4.16) Mij =0, 1 < i < j < n - 1. 
Restriction of the system on invariant manifold is given with: 

n-1 

M in = -J2^Yj M jn-XnTi, t = l,...,n-l, 
j=i 

n-1 

(4.17) f i = -(Ji + J„)r„M m - nVTj, i = 1, . . . ,n - 1, 

n-1 

Tn - (Ji + J„)^L,M Jn , 

i=i 

w/iere ft n = M J n + J n M = -J ln E^M^i A E t + M U E, A E n ). 

Remark 4.1. It can be proved that equations (4.10) have an invariant relation (4.16) if 
and only if pr B oA o pr B is proportional to the identity operator on d [36] . 

Applying a general construction of compatible Poisson structures related to the sym- 
metric pair decomposition of semi-simple Lie algebras for the symmetric pair (so(n + 
1), so(n)) (sec [66, 70, 10]), one gets the second Poisson structure on e(n): 

(4.18) {r,,M,„} 2 = 5ij, {TuTjh = -Mij. 

Let us denote x„ = a = Ji + J n = 1, 6 = 2J\. 

Proposition 4.2. (i) The Casimirs of the structure (4.18) are functions (4.14) together 
with I = \{Ml n + M| n + • • • + M2_ 1; J + r„. 

(ii) The Lagrange top is Hamiltonian with respect to the Poisson structure (4.18) and 
the Hamiltonian function: 

1 n — 1 

h a , 2 = -(r,r) + (b- 1) ^iMijM jn . 
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In this setting, the system is integrable in the commutative sense by means of integrals 
(4.13). The invariant tori are of dimension n — 2. 

Therefore, the system (4.10), (4.3), (4.15) satisfies axioms (Al), (A2'), (HP) and 
(BP). 

The restricted Hess-Appel'rot system (4.17) admits the Lax representation (4.12) 
with M and f2 related by (4.3), (4.15). But, only 3 integrals from the family (4.13) are 
independent when invariant relations (4.16) are satisfied. Namely, the spectral curve is 
then given by [36]: 

(4.19) p(X, n) = det(L(A) - /Ad) = (-^)"~ 3 (n* + (J?PW + 0(A)) = 0, 
P(A)^F 2 + ^A 2 F 1 + A 4 (^) 2 , 0(A) = A 2 P 3 , 

where Fi, F 2 , F 3 are integrals 

n—l n 

Fi=H=^Y, M ?n + XnT n , F 2 = ^r 2 = l, 

i=l i=l 

(4.20) f 3 = Yl (M m r, -M,„r 4 ) 2 , 

l<i<j<n-l 

that correspond to the integrals (2.2) of the classical problem. 

On the invariant submanifold (4.16), integrals Fi,F 2 ,F 3 are not enough to give an 
integrability of the Lagrange top. But, when conditions (4.16) are satisfied, the system 
(4.10), (4.11) has additional integrals 

Fij = M in Tj — Mj n Ti, l<i<j <n-l 

that implies that the Lagrange top on invariant submanifold (4.16) can be solved by 
quadratures: (4.16) is almost everywhere foliated on invariant 2-dimensional tori. 

Integrals F^ are not integrals of Hess-Appel'rot system (4.17). Nevertheless, the 
system (4.17) has additional invariant relations 

(4.21) F^ = M in Tj - M jn Ti = 0, 1 < i < j < n - 1, 

which are equivalent to condition F 3 = 0. Relations (4.21) show collinearity of 
vectors (M ln ,M 2n , . . . ,M„_i,„) and (ri,r 2 , . . . ,L„-i)- 

4.3. Mishchenko Fomenko flows. Let G be a compact Lie group, g its Lie al- 
gebra, (•, •) a Adc-invariant scalar product on q. Let a € Q be an arbitrary element 
Q a = {v G g, [a, rf\ = 0} be the isotropy algebra and G a be the adjoint isotropy group of 
the element a. 

Let g — Q a © f be the orthogonal decomposition. Consider the linear operator (so 
called sectional operator [70]) A a ^^ ■ — > 0, defined by 

(4.22) A aAC (0 = ad" 1 o ad b o pr B (£) + C(pr fla 0, 

where b belongs to the center of a , pr and pr fla are the orthogonal (with respect to (•,•)) 
projections to and a , respectively and C : Q a — > a is a positive definite, symmetric 
operator such that the quadratic form (£, C{£)) is Adc Q -invariant: 

(4.23) [£,C(0]=0, Ce 0a . 

We can always find b and C such that A a ^^c is positive definite. If a is regular, i.e., 
G a is commutative, then the condition (4.23) is always satisfied. 
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Identify g* with g by means of the scalar product (•, •) and consider the left-triviali- 
zation: T*G =i G x g = {(g,£,)}- Then the quadratic form 

(4-24) H aAC = ^(A aAC (0,0 

can be regarded as the Hamiltonian of a left-invariant Riemannian metric K a .b.c on G. 
After left G-reduction, the equations of the geodesic flow take the form of Euler equations 

(4-25) i = % Vfr oA c(0] = [£, A aAC (Ol 

that are Hamiltonian with respect to the Lie-Poisson bracket 

(4.26) {F 1 ,F 2 } 1 = -(^[VF 1 ,VF 2 }}. 

The Casimir functions of the Lie-Poisson brackets are invariant polynomials. 

Let S be the set of linear functions on g a . Let pi, . . . ,p ra nkG be the basic invariant 
polynomials of the Lie algebra g. Mishchcnko and Fomcnko proved that the polynomials 
obtained by shifting of argument of invariants: 

(4.27) C: =Pi(£ + \a), i = 1, . . . , rankG, AeM 

are commuting integrals of the system (4.25) and {C,S}i = [52]. The algebra C + S 
is a complete algebra integrals of the system (4.25) [52, 10, 70]. We refer to systems 
(4.25) as Mishchenko-Fomenko flows. 

Moreover, the system is bi-Hamiltonian (see [51, 70]). The second Poisson structure 
is linear: 

(4.28) {F 1 ,F 2 } 2 = -{a, [VFi, VF 2 ]>, 

with a set of Casimirs S. The integrals Pi t \(£), i = 1, . . . ,rankG are Casimir functions 
of the bracket {•, -}i + A{-, -} 2 . 

Now, let us perturb the Hamiltonian (4.24) as follows: 

(4-29) H aAC >,D = \(A aAC >(Z),0 + (D(pr,0.0. 

where D : D — > g a and A aA c' is given by (4.22) such that H aA c',D is positive definite. 
Here C not need to satisfy (4.23). 

Proposition 4.3. The perturbed system 

(4.30) C = K,vzr o , 6> c,z>(0] 

has an invariant manifold 

(4-31) Pr Sa (£) = ° 

and satisfies axioms (Al), (A2), (HP), (BP) of the systems of the Hess-Appelrot type. 



5. Partial reductions 

5.1. Classical Hess-Appel'rot system and spherical pendulum. Let us return 
to the classical problem (2.1), (2.5). Consider the motion of the rigid body in the space 
reference frame. Let R <G 5*0(3) be the matrix that maps the moving reference frame to 
the fixed one. Following Arnol'd's notation [6], denote m = i?M,w = i?f2,7 = RT. The 
position of mass centre is z e 3 , where e 3 = i?E 3 = i?(0, 0, 1) T . 
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Beside geometric interpretation of the conditions (2.3), Zhukovski also noticed (see 
[73]) that the motion of e 3 is described by the spherical pendulum equation. Indeed, the 
Euler-Poisson equations for the variables (m, e 3 ), in the space frame, read 

(5.1) m = z 7xe 3 , 

(5.2) e' 3 = u> x e 3 = pr e x uj x e 3 . 

On the invariant set M 3 = we have pr E J- fl = (0,1,0,2,0) = (J\Mi, J\Mi, 0) = 
JiM. The relation in the space frame gives 

(5.3) pr e j_ w = Jim. 

By differentiation of (5.2), using (5.1) and (5.3), we get 

e 3 = Jiz (-f x e 3 ) x e 3 + jfm x (m x e 3 ) = -Jiz j + e 3 ( Jiz (e 3 , 7) - (e' 3 ,e' 3 )), 

describing the Euler-Lagrangc equations with multiplier A = Ji^o(e 3 ,7) — (e' 3 ,e' 3 ) on 
the sphere (e 3 ,e 3 ) = 1. This is the pendulum system with Lagrangian 

l(e 3 ,e 3 ) = ^j-( e 3,e 3 ) -z (7,e 3 ). 

In particular, the motion of e 3 can be found by elliptic quadratures. 

5.2. Invariant relations and reductions. Note that, considered on the whole 
phase space T*SO(3) of the rigid body motion, the function M3 is the momentum map- 
ping of the right S'0(2)-action-rotations of the body around the line directed to the center 
of the mass. By the analogy with the reduction of the system to the spherical pendu- 
lum, in this subsection we study reductions of the Hamiltonian flows that satisfy axioms 
(Al), (A2), (HP) of the systems of the Hess-Appelrot type restricted to their invariant 
submanifolds. Apparently, the lowering of order in Hamiltonian systems having invariant 
relations was firstly studied by Levi-Civita (e.g., see [41, ch. X]). 

Let G be a connected Lie group with a free proper Hamiltonian action on a symplectic 
manifold (M,u) with the momentum map (9.9). Assume that is a regular value of 
Let £1, . . . , £ p be the base of g. Then the zero level set of the momentum mapping (9.9) 
is given by the equations 

(5.4) M : & = = i = l,...,p. 

Let (N n ,ujo) be the symplectic reduced space and let 7r : M — > iVo = M /G be the 
canonical projection (see subsection 9.4). Consider the Hamiltonian equations: 

(5.5) x = X h (x). 

Theorem 5.1 ([36, 37]). (i) Suppose that the restriction of h to (5.4) is a G-invariant 
function. Then Mq is an invariant manifold of the Hamiltonian system (5.5) and Xh\M a 
projects to the Hamiltonian vector field Xh : dno(Xh)\x = Xh [wx); where h is the 
induced function on N defined by 

(5.6) h\ M „ = TToh v = h o tt . 

(ii) The inverse statement also holds: if (5.4) is an invariant submanifold of the 
Hamiltonian system (5.5), then the restriction of h to Mq is a G-invariant function and 
Xh\M projects to the Hamiltonian vector field Xh on N , where h is defined by (5.6). 

In both cases, the Hamiltonian vector field Xh is not assumed to be G-invariant on 
M. Moreover Xh\M may not be G-invariant as well. It is invariant modulo the kernel 
of d7r , which is sufficient the tools of symplectic reduction are still applicable. 
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We shall refer to the passing from x = Xh\M to 

(5.7) y = X ho 

as a partial reduction. 

The natural way to obtain a Hamiltonian h such that corresponding flow have in- 
variant relations (5.4) is to perturb a G-invariant Hamiltonian h\ G Cq{M) 

p 

(5.8) h = h K + Y J h i (j>i, 

i=i 

where hi are arbitrary smooth functions on M. Then the set of p+1 functions h, <j>\ , . . . , (f>k 
satisfy axioms (Al), (A2'), (HP) of the systems of the Hess-Appelrot type, i.e., the Hamil- 
tonian system (5.5) is restrectively integrable and of the form of Hamiltonian perturbation 
(see [21]). 

The function h\ is G-invariant and we can perform the usual symplectic reduction 
to the Hamiltonian flow on the symplectic reduced space No. Since the Hamiltonians h 
and Ha coincide on M , the reduced flow is the same as the partially reduced flow (5.7). 
However, Hamiltonian vector fields Xh\M an d A/j a |m ar e different. 

The partial reduction can be seen as a special case of the symplectic reductions 
studied in [8, 43] (see also [42], Theorem 14.6, Ch. III). 

Theorem 5.2 ([42]). Let (M,oj,h) be a Hamiltonian system and let M C M be an 
invariant submanifold of the vector field Xh upon which the symplectic form w induces a 
2- form u>m of constant rank. Assume that there is a surjective submersion with connected 
fibres 7To : Mq — ► Nq onto another symplectic manifold (Nq,u!q), which satisfies tTqLUo = 
u)\m - Then there exist a unique reduced Hamiltonian function ho on N such that H\m = 
ho o7T and dir (X h )\ x = X ho \^ o{x) . 

5.3. Zhukovskii property. An immediate corollary of theorem 5.1 is 

Corollary 5.1. Suppose that the partially reduced system (5.7) is completely integrable 
in the non-commutative sense and No is almost everywhere foliated on r-dimensional 
isotropic invariant manifolds, level sets of integrals ff , i = 1, . . . ,dim7Vo — r. Then M 
is almost everywhere foliated on (r + dim G) -invariant isotropic manifolds 

(5.9) M c = {fi = 7r*.f° = Ci | * = 1, ... .dimJVo - r} 
of the system (5.5). 

In the case of the classical Hess-Appel'rot system we have M = T*SO(3), G = 50(2) 
(rotations of the body around the vector E 3 ) and the reduced system is the spherical 
pendulum. The reduced phase phase N = T*S 2 is foliated on two-dimensional tori, while 
the invariant manifold Mo C T*S0(3) defined by M3 = is foliated on 3-dimensional 
invariant Lagrangian tori. 

Definition 5.1. We shall say that the Hamiltonian system (5.5) has the Zhukovskii 
property if it has invariant relations of the form (5.4) and that the reduced system (5.7) 
is completely integrable. 

Let us make a terminological note: in [36, 37] for the Zhukovskii property has been 
used a different name, partial integrability (or geometrical Hess-Appol'rot conditions for 
natural mechanical systems). However, there is an another notion named partial in- 
tegrability which has been introduced in a study based on the Poincare-Lyapounov- 
Nckhoroshev theorem [60, 32, 61] (see also [33, 45]). There, a Hamiltonian system 
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(5.5) restricted to an invariant submanifold N C M of lower dimension is completely 
integrable, i.e., N is filled with periodic or quasi-periodic trajectories of (5.5). 

The following example will illustrate the difference between the usual and partial 
reduction of Hamiltonian flows. 

Example 5.1. Let G be a torus T p and let the reduced flow be completely integrable 
in the commutative sense by integrals /°, . . . , m = \ dim N . Consider the regu- 
lar compact connected component level set of /°, . . . , By Liouville's theorem, it is 
diffcomorphic to a m-dimcnsional torus T m with quasi-periodic flow of (5.7). Thus the 
compact connected component M c = Tr^ 1 (T m ) of (5.9) is a torus bundle over T m : 

T p — > M c 

TO 

Suppose that fi,. ■ ■ ,f m can be extended to commuting T p -invariant functions in 
some T p -invariant neighborhood V of M c - Then, within V, M c is given by the equations 

fx = ci, . . . , f m = c TO , fa = 0, . . . , <f> p = 0. 

From the Noether theorem the functions ^ commute with all T™-invariant functions 
on M and the following commuting relations hold on V: 

{fajb} = {fa,4>i} = {<!>i,<j>j} = 0, 
a,b=l,...,m, i,j = l,...,p. 

Now, as in the case of commutative integrability of Hamiltonian systems M c is a 
Lagrangian torus with tangent space spanned by Xf a , X^ t , i.e., the bundle (5.10) is 
trivial. However, in general, the flow of Xh over the torus M. c is not quasi-periodic: the 
vector field Xh does not commute with vector fields 

X fl , . . . , X/ m , , . . . ,X^ v 

although Poisson brackets {h, /„}, {h, fa} vanish on M c - So, in general, we can not apply 
the Lie theorem [39] to solve the system by quadratures. 

Note that, if h is G-invariant, the complete integrability of the reduced and original 
system are closely related [74, 37]. In particular, if as above G is torus, then the re- 
construction problem is easily solvable by quadratures (e.g., see [48]). By contrary, in 
the case of partial reductions, although the reduced motion y(t) is found by quadratures, 
the reconstruction problem in determining x(t) (no(x(t)) = y(t), x(to) G ttq (y(to))) in 
general, leads to non-autonomous equations. 

5.4. Reductions of additional symmetries. Suppose that an additional free 
Hamiltonian action of a connected compact Lie group K is given. Let 

* : M -> fi* 

be the corresponding momentum mapping and let 

a : M -> M/K 

be the natural projection. In what follows, by capital letters we shall denote the functions 
on P = M/K and with small letters corresponding X-invariant functions on M: 

a* :C°°(P) -^C^(M), F h-> / = F o a. 



(5.10) 
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Since the action of A is Hamiltonian, the Poisson bracket of two if -invariant func- 
tions is if -invariant as well, so manifold P = Mj A carries the induced Poisson structure 
{;-} K defined by a*{F 1 ,F 2 } K = {a*F u a*F 2 }. 

The symplcctic leaves in (P,{-,-} K ) are of the form 1 i'~ 1 (O v )/K, where O v is a 
coadjoint orbit of r/ 6 t*. The Casimir functions I\, . . . , I r of (P, {•, - } K ) are: 



i j = (<r*r 1 (p j on 



, r = rank A, 



where p\,. . . ,p r is the base of homogeneous invariants on f,r = rank if. 

Further, suppose that the actions of G and A commute, that is we have ipj} = 0, 
i = 1,- --,P, j = 1,- --,q, where 



^ = (*,&). 



1. 



,q. 



and Ci, • ■ • , C q is a base of t. Thus </>i are if- invariant and we have induced G-action on 



P given by the Hamiltonian functions $j, <&j = (er*)~ 



l,...,f>. 



Suppose the functions h^,h\, . . . ,h p in the definition of the perturbed Hamiltonian 
(5.8) are A-invariant functions. Then we can reduce the system (5.5) to the Poisson 
manifold P as well: 

(5.11) F = {F,H} K , FeC°°(P). 

Here the reduced Hamiltonian H is of the form 



(5.12) 



where H\ is G-invariant, i.e., {iiAj'l'i}^ 
relations 



0. As a result we obtain reduced invariant 



(5.13) 



*i = 0, i = l, 



defining the invariant manifold Po = M /K of the reduced flow (5.11). 

The functions $i are integrals of the non-perturbed flow so H satisfy axioms (Al), 
(A2') and (HP) of systems of Hess-Appel'rot type. 

On the other hand, we have the induced Hamiltonian A-action on (A ,o>o) with the 
momentum mapping 'J'o satisfying ^>\m = $o OI o (see diagram (5.14) below). Since 
h is A-invariant, the reduced Hamiltonian h on N n is also A-invariant. Therefore the 
momentum mapping *S? is conserved along the flow of the partially reduced system (5.7). 

r r e* 



(5.14) 0*^-M D Mo = $ _1 (0) ^ N = M /G 

1 1 

P = M/K D P a = Mq/K 

Remark 5.1. Suppose we have additional integrals of the reduced system (5.7) implying 
Zhukovskii property of (5.5) with respect to the G-action. Then Mq is almost everywhere 
foliated on invariant isotropic manifolds M c (sec Corollary 5.1). Hence we obtain invari- 
ant foliation of Po = Mo /A on manifolds of the form V c — (A • M. c )/K. It is clear that 
V c are isotropic submanifolds of the appropriate symplectic leaves in (P, {•, - } K ). 
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5.5. Natural mechanical systems. Let (Q, k, v) be a natural mechanical system, 
where the metric k is also regarded as a mapping k : TQ — > T*Q. Let G be a connected 
Lie group acting freely and properly on Q and p : Q — > Q/G be the canonical projection. 

For Lagrangian systems, it is convenient to work with tangent bundle reductions. 
Let V q = {£ q £ £ fl} be the tangent space to the fiber G ■ q (vertical space at q) and 
V = U q V q be the vertical distribution. Consider the horizontal distribution H orthogonal 
to V with respect to the metric k: 

(5.15) H = {X q e T q Q I ( Kq (X q ), = 0, £ e fl, g € Q} 

Let $ be the cotangent bundle momentum mapping (9.10). Since W = k _1 ($ _1 (0)), 
the horizontal distribution is invariant with respect to the "twisted" G-action 

(5.16) go(q,X) = (g-q, K - q o(dg- 1 )*o Kq (X)), XeT q Q, 

that is the pull-back of canonical symplectic G-action on T*Q via metric k. With the 
above notation we get [36] 

Theorem 5.3. (i) (Partial Noether theorem) The horizontal distribution (5.15) is an 
invariant submanifold of the Euler- Lagrange equations (9.4) of the natural mechanical 
system (Q, k, v) if and only if the potential v and the restriction k-h of the metric k to H. 
are G-invariant with respect to the action (5.16). 

(ii) (Partial Lagrange-Routh reduction) If TL is an invariant submanifold of the 
system (Q,n,v), then the trajectories q(t) with velocities q(t) that belong to H project to 
the trajectories Tt(q(t) of the natural mechanical system (Q/G,K,V) with the potential 
V(ir(q)) — v(q) and the metric K obtained from k-h via identification TL/G = T(Q/G). 

Note that when k is G-invariant, the twisted G-action (5.16) coincides with usual 
G-action: g ■ (q,X) = (g ■ q,dg(X)) and the induced metric K is the submersion metric. 

6. Systems of Hess Appel'rot type and Zhukovskii property 

In the sequel we shall study the relationship between the Zhukovskii property and 
systems of Hess-Appel'rot type on the Poisson manifold (P,{-,-} K ). In particular, we 
shall establish the dimension of invariant manifolds stated in Remark 5.1. 

We say that some property holds for a generic point x of the manifold M if the 
property holds on an open dense set U C M. 

6.1. Hamiltonian perturbation. Consider the system (5.11) with Hamiltonian 
function H of the form (5.12), where Hi are arbitrary smooth functions on P and Ha G 
Cq(P) is a G-invariant Hamiltonian. 

Lemma 6.1. If F is a G-invariant integral of the system 

(6.1) F = {F,H A } K , FeC°°(P), 

then it is also a integral of the system (5.11) with perturbed Hamiltonian function (5.12), 
restricted to the invariant manifold (5.13). 

Proof. Directly, from {$i,F } K = and {H Al F} K = we get 
{H,F} K \ Po = {H A ,F} K \ Po +j2{H^,F} K \ Po 

i=l 

= H ^ F j} K \?o + E F ^\Po = °- 

i=l i=l 

□ 
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Therefore, if the non-perturbed flow (6.1) is integrable, for a generic perturbation 
(5.8), only G-invariant integrals remain to be integrals of the perturbed flow restricted 
to P . 

In what follows, we suppose that at a generic point z of Pq, functions I\, . . . ,I r are 
independent, or equivalently, the symplectic leaf through z is regular. 

Let F\, . . . , F p be the Poisson-commuting integrals of the non-perturbed system (6.1), 
mutually independent and independent of functions I\, . . . , I r , $i, . . . , $ p at a generic 
point z e P . Let 

s = dimG z = dim((rf/i, . . . , dl r ) n (d<&i, . . . , d& p ))\ z , 

at a generic z e P (G z C G is the isotropy group of z with respect to the induced 
G-action) and let 

(6.2) T ={*i,... ,%,h, ... ,I r ,F u ... ,F p } c C°°(P), 

(6.3) T a = {<, ■ • • , <, f°, • • • , f p } C G°°(A ). 

Here if)® = (^o, Q) are component of the momentum mapping ^Sf and f are obtained 
from Fi by the composition 

(6.4) C%(P) Gg^(M) G£(M ) G^(iV ), 
where i : M ^ M is the inclusion. 

Theorem 6.1. (i) 

(6.5) 2p < dimP + 2s - 2dimG - rank if. 

(ii) If T is a complete set at a generic point z G P then G is commutative and the 
partially reduced flow (5.7) is completely integrable by means of integrals Tq. 

(hi) Contrary, suppose Tq is a complete set of functions on Nq. Then the manifold 
Pq is almost everywhere foliated on i (dim P — rank K) -dimensional level sets of integrals 
T , both of the non-perturbed (6.1) and the perturbed system (5.11). This is enough for 
integrability of the non-perturbed system (6.1) in the case when G is a commutative Lie 
group. 

PROOF. I Consider Theorem 9.1 where we put (M, u, G, $) to be equal to (No, uj a , K, *o)-| 
We have (see [11]): 

dind y* C°°(F) = dind C^(N ) = 

(6.6) = dind (**C°°(r) + C£(N )) = dimi^ - dimK y 
ddim G£? (N ) = dim N a + dim K y - dim K 

for a generic y £ N , p = ^fo(y). 

As above, let pi, . . . ,p r be the base of homogeneous invariants on F, r = rank if. 
Then functions 

i9= Pi o^ , i = l,...,r 
belong to $>*C°°(t*) n C^(N ), and therefore belong to the centers in both algebras 
**G°°(r) and C%(N ). Among them there are exactly dind C^(N ) independent ones. 
Under the composition (6.4), corresponds to the Casimir function Ij of the bracket 

{■,■}*■ 

Since F\ , . . . , F p are mutually independent and independent of the functions I\,...,I r 
and $i, . . . , $ p at a generic point z e F , it follows that f, ■ ■ ■ , f are also mutually 
independent and independent of the functions i\, . . . , i® and 

(6.7) ddim {*°, ••• ,i° r ,f,..., fp} = dind G^(iV ) + p. 
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The set {i\, ■ ■ ■ , i®, /i , ■ • ■ , fp} is a commutative subset of C^(N ), so 

(6.8) ddim {*°, ••• ,^,/i ,..., f p } < i(ddim G£(jV ) + dind C%(N )) 

The condition that a symplectic leaf through z £ Po is regular implies that the 
coadjoint orbit through p = ^o{y) is regular, where y = no(x), z — a(x), x G Mq. 
Therefore, from (6.6) we get 

ddim C% (N Q ) + dind C% (N ) = dim N + rank K - dim K 

(6.9) = dimP-2dimG + rank/s: 
Also s — dimG^ = dimK y , for a generic y = tt (x), z = a(x), x e M . Thus 

(6.10) dind C%(N ) = rankPJ - s . 

Finally, combining (6.7), (6.8), (6.9) and (6.10) we obtain inequality (6.5). 

II By using the proof of item (i), we get 

ddim T = dim G + rank K + p — s, 
dind T = rank G + rank K + p — s, 

ddim T + dind T = dim G + rank G + 2 rank if + 2p - 2s 

(6.11) < dim G + rank G + 2 rank if - 2s + dim P + 2s - 2 dim G - rank K 

(6.12) = dimP + rankG - dimG + rank AT < dim P + rank K. 

If (6.2) is a complete set at a generic point z e P then we have equalities both in 
(6.11) and (6.12). Equality in (6.12) implies that G is a commutative group. Equality in 
(6.11) means that we have equality in (6.5), that is, according (6.8), {i\, ■ ■ ■ , i®, /i , • • • , fp} 
is a complete commutative subset in C^(N ). 

Then, from Theorem 9.1 we get that {ip®, . . . , ipq,i1, ■ ■ ■ , i®, f\ , ■ ■ ■ , fp} is a complete 
set on No- Since i° are polynomial functions in ip®, the set (6.3) will be also a complete 
set of integrals of the partially reduced system (5.7). 

III Suppose (6.3) is a complete set. Then {i®, ■ ■ ■ , i®, fi , ■ ■ ■ , fp} is a complete 
commutative subset in C^(N ) and we have equality in (6.5). Therefore ddim T = 
dimG + rankif + ±(dimP + 2s - 2dimG - rankif) - s = ±(dimP + rankif). The 
dimension of invariant level-sets given by T is dim P — ddim T = \ (dim P — rank K) . □ 

6.2. Lifting of bi Poisson structure. As above, consider the system (5.11) with 
Hamiltonian function H of the form (5.12), where Hi are arbitrary smooth functions 
on P and Pa S Cq(P). Suppose in addition to the Poisson structure {•, - }i = {-,-} K i 
that the non-perturbed system (6.1) is Hamiltonian with respect to the another Poisson 
structure {•, -}2 which is compatible with the first one. Also, we suppose that functions 
$i are Casimir functions of the second bracket: 

(6.13) {$ 4 ,P} 2 =0, i = l,...,p, PeG°°(P). 

Thus, the Hamiltonian flow (5.11) satisfies axioms (Al), (A2'), (HP) and (BP) of 
systems of Hess-Appel'rot type. 
Let 

n = {{•, -}x lM = Ai{-, '}i + -h I Ai, A 2 e R, A? + \% + 0}, 
be the corresponding pencil of compatible Poisson structures. 

Theorem 6.2. The pencil of compatible Poisson structures U on P induces the pencil 
of compatible Poisson structures Hq within the algebra Cj-p(iVo) of K -invariant functions 
on No 
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PROOF. The Poisson brackets {-,-}i and {-,-} 2 are compatible. That is why it is 
enough to prove that both Poisson brackets are correctly defined within the algebra of 
functions C%{N ). 

The G and K actions on (M,ui) is Hamiltonian. Thus {•, -}i is well defined in all of 
the algebras: Cg(M) = C°°(P), C°°(M), Cg(M), C°°(N Q ) and C%(N ). 
The function F G G°° (P) is G-invariant if and only if 

(6.14) {F,$i}i=0, i=l,..., p. 

The bracket {■, -} 2 induces the Poisson structure within Cq(P) if the bracket {Fi, F 2 } 2 | 
of two G-invariant functions is again a G-invariant function. 

Suppose Fi and F 2 are G-invariant. Let us write the Jacobi identity for the bracket 
{•7 }i + {•> and functions Fi,F 2 , <&,: 

{{F 1 ,F 2 } 1 + {F 1 ,F 2 } 2 ,$ i } 1 + {{F 1 ,F2} 1 + {F 1 ,F 2 } 2 ,$ i }2 + 
Fiji + Fi} 2 , F 2 } 1 + {{$,, Fiji + Fi} 2) F 2 } 2 + 
{{F 2 , + {F 2 , $J 2 , Fi}! + {{F 2 , $j! + {F 2 , $i} 2 , Fi} 2 = 0. 

Taking into account identities (6.13) and (6.14) for $ i; Fi, F 2 , from the Jacobi identity 
we get 

(6.15) {{Fi,F 2 }i,$Ji + {{Fi,F 2 } 2 ,$Ji = 0. 

Since Fi and F 2 are G-invariant, {Fi, F 2 }i is also G-invariant, so the first term is equal 
to zero. Hence, the second term in (6.15) equals zero as well and the bracket {•, -} 2 
induces the Poisson structure within Cg>(P) = Cg G {M). 

The functions <£>i are Casimir functions for the bracket {-,-} 2 . Whence, after the 
restriction, we have well defined bracket {•, -} 2 on the submanifold (5.13). From the 
above considerations it follows that the Poisson bracket {•, -} 2 is well defined within the 
algebra of functions G^(F ) = C^ K {M ) C£(N ). ' ' □ 

The pencil LTo can be used as a tool in proving complete intcgrability of the partially 
reduced system (5.7). For example, suppose that all structures in II, except possible 
{■, -} 2 , have the maximal rank equal to dim P — rankF, and that their Casimir functions 
are globally defined on P. Then the union C of Casimir functions of all Poisson structures 
from IT non-proportional to {•, -} 2 is a commutative set with respect to the all Poisson 
structures from II. Moreover, let S be the set of Casimir functions of the bracket {•, -} 2 . 
Then {C 7 S}i = and, since {$1, . . . , & p } C S, the functions in C are G-invariant. If 
the pencil II satisfies conditions of the form (9.11) (for more details see [10]), the set 
of functions T = C + S will be a complete set of functions. However, as we have seen 
in Theorem 6.1, even if T is complete, it not need to be complete at the points of the 
invariant set Pq. 

Nevertheless, suppose that the Casimir functions C induce the Casimir functions 
Co C Ck(No) of the brackets non-proportional to {•, -} 2 . Then, from Theorems 9.3 and 
9.1 we get 

Proposition 6.1. The set C * s a complete commutative subset of C^(N ) if and only 
if the corank of the coplexified bivectors 

^ : (ann(F • y)f x (ann(F • y)f - C 

is equal to dind C^(Nq) for all tt G Hq at the generic point y G Nq. In this case the 
partially reduced system (5.7) is integrable by means of integrals Tq = Co + {ipi, ■ ■ ■ , ipq}- 
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Here axm{K ■ y) C T*N is the annihilator of the tangent space to the orbit K ■ y of 

y- 

The condition stated in Proposition 6.1 is verified in proving the complete commu- 
tative integrability of geodesic flows of normal metrics on adjoint orbits in [12, 58] (see 
Remark 7.2 given below). 

7. Zhukovskii property. Examples 

7.1. Magnetic flows on adjoint orbits. In order to describe partial reduction 
of the Hess-Appelrot rigid body system (4.1), (4.3), (4.4), (4.7) and the geodesic flows 
(4.30) we shall need a description of certain natural mechanical systems on adjoint orbits 
recently studied, in the presence of the additional magnetic force, in [13, 14]. Note that 
integrable magnetic flows on homogeneous spaces are also given in [66, 28, 46]. 

Let G be a compact connected Lie group with the Lie algebra g and invariant scalar 
product (•, •). Consider the G-adjoint orbit 0(a) C fl. The tangent space at x = Ad g (a) 
is simply the orthogonal complement to g x = {£ G g [x,£] = 0}. The cotangent bundle 
T*0(a) can be represented as a submanifold of g x g: 

T*0(a) = {(x,p) | x = Ad g (a),p G g^}, 

with the paring between p G T*0{a) = g^ and 77 G T x O(a) given by p(rj) = (p,r]). 
Then the canonical symplectic form u> on T*0(a) can be seen as a restriction of the 
canonical linear symplectic form of the ambient space g x g: Y^i=i 3 dpi A dxi, where pi, 
Xj are coordinates of p and x with respect to some base of g. Let il be the standard 
Kirillov-Kostant symplectic form on 0(a). 

The canonical G-action g ■ (x,p) = (Ad g x, Ad g p) on the magnetic cotangent bundle 
(T*0(a),ui + ep*fl) is Hamiltonian with the momentum mapping (see [28, 13]) 

(7.1) m : T*0(a) -> g* S g, m{x,p) = [x,p] + ex. 

A natural generalization of the magnetic spherical pendulum to the orbit 0(a) is the 
mechanical system with the kinetic energy given by the normal metric and the potential 
function V(x) = (c, x), i.e., with the Hamiltonian 

h c {x,p) = ^([x,p], [x,p]) + (c,x). 

Another natural class of systems are the magnetic geodesic flows of the G-invariant 
metrics K a ^ defined by the Hamiltonian function 

h a ,b( x >P) = ^([bxP], [x,p]) = -^(&d x a,d bx p,p), 

where b belongs to the center of g a and b x = Ad g (b) for x = Ad g (a) [13]. For compact 
groups, we can take b such that K a ^ is positive definite. If b = a we get the Hamiltonian 
of the normal metric. 

The equations of the magnetic pendulum, in redundant variables (x,p), are given by 
(see [14]) 

x = [x,\p,x]], 

P= [P, [P, x]] + e[x,p] - c + pr 0x c, 
while the magnetic geodesic flow read (see [13]) 
x = -ad K ad bx p = [[6 x ,p],x], 

(7.3) 

p = -ad,: 1 ^, [x, [b x ,p]}] +pr [[b x , p},p] + e[b x ,p\. 
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Let 0(a) be an arbitrary orbit. 

Theorem 7.1 ([13, 14]). The magnetic pendulum system, for a regular element c£g, 
and the geodesic flows of the metrics K a ^ on (T* O(a) , u> + ep*£l) , described by equations 
(7.2) and (7.3), respectively, are completely integrable by means of polynomial integrals. 

7.2. Partial reduction of rigid body systems. Consider the construction given 
is Section 4 for the case when the symplectic manifold M is the phase space of the n- 
dimensional rigid body motion about a fixed point T*SO(n). As usual, we use the left 
trivialization T*SO(n) =i SO(n) x so(n)(g, M), where so{n) and so(n)* are identified 
by the use of invariant scalar product (X, Y) = — ^tr(XY). 

We follow Ratiu's generalization of the heavy rigid body motion [64] . Let 

7 = 7i2^i AE 2 + 734^3 A E 4 + ■ ■ ■ + lk ^ k E k ^ A E k , k = 2{n/2] = 2 rankSO(n), 
where 7i,i+i are mutually different. Also, let 

X = X12-&L AE 2 , n>4, 
X = X12-B1 A E 2 + X34-E3 A E 4 , n = 4, 

where X12 7^ X34> X12 7^ 0. Then the adjoint isotropy groups of 7 and x and corresponding 
isotropy Lie algebras read, respectively (e.g., see [9]) 

k times 

K = SO{n) y = 50(2) x •• • x SO(2), t = {E x A E 2 , ■ ■ ■ , £ fe -i A E k ), 

G = SO(n) x = SO(2) x SO(n - 2), g = (£4 A E 2 , Ei A Ej | 3 < i < j < n). 

Note that if is a maximal torus in SO(n), so the adjoint orbit through 7 is the flag 
manifold 0(7) = SO(n)/K, while the adjoint orbit 0(x) = SO(n)/G = SO(n)/SO(2) x 
SO (n — 2) is Grassmannian variety Gr + (n, 2) of oriented 2-dimensional planes through 
the origin in M™. 

Consider the natural left action of K and right action of G on T*SO(n). The corre- 
sponding momentum mapping ^> and <f>, in the left trivialization are given by 

*(g, M) = pr { (Ad 5 M), <%, M) = pr g (M) 

The Hamiltonian of the n-dimensional Hess-Appelrot rigid body system is the Hamil- 
tonian perturbation of the Lagrange system (or Lagrange bi-top system for n = 4): 

(7.4) h(g,M) - l(M,fi> + (Ad fl -i 7 ,X>, 

where il = JM + MJ and J is given by (4.7). The fixed element 7 play the role of 
the horizontal vector (direction of the gravitational force) as seen in the space reference 
frame. 

Lemma 7.1. The Hamiltonian (7.4) is left K-invariant and right G-invariant on the 
zero level set of the momentum mapping $; 

(7.5) (T*SO(n)) Q : 4>i 2 {M,g) = M 12 = 0, ^(M,g) - My = 0, 3 < * < j < n. 

We have P = (T*SO(n))/K = so(n) x 0(j)(M, L). The X-reduced flow is described 
by equations (4.1), (4.3), (4.7) (or (4.4) for n = 4), where one have to fix values of 
invariants in T in order that T belongs to the adjoint orbit 0(7) 

On the second hand, the symplectic reduced space No = (T*SO(n))o/G is symplec- 
tomorphic to the cotangent bundle of the adjoint orbit 0(x) = Gr + (n, 2). 
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Let x = Ad g x £ 0(\) (in 3-dimensional case, x represents the position of the mass 
center in the space coordinates). We can rewrite the Hamiltonian (7.4) in the form 

h(M, 9)= Jl+ 2 Jn (M, M) + 12 H 12 + fcjHii + ^x). 

3<i<j<n 

Therefore 

h(M,g)\ (T , so{n))o = •I±±^(M,M) + < 7 ,x). 

After reduction to T*0(x), the bi-invariant kinetic energy term goes to the kinetic 
energy of the normal metric multiplied by (Ji + J n ) (e.g, see [9]). Thus, likewise in the 
3-dimensional case, the partially reduced flow is the pendulum system on O(x)- 

x = (Ji + J n )[[x,p],x], 

P = (Jl + J n )[[x,p],p] - 7 + P r so(«)* T' 

with Hamiltonian h (x,p) = Jl + J " ([x,p], [x,p]) + (x, 7) (we follow the notation of the 
previous section). It follows from Theorem 7.1 that the reduced system is completely 
integrablc. Hence the system satisfies Zhukovskii property and we get the following 
qualitative behavior of the system, 

Theorem 7.2. The partial reduction of the Hess-Appel'rot rigid body problem defined by 
the Hamiltonian (7.4) is completely integrable pendulum type system (7.6) on the oriented 
Grassmannian variety Gr + (n, 2). The invariant manifold (7.5) is almost everywhere 
foliated by invariant dim SO(n)- dimensional Lagrangian invariant manifolds that project 
to the 2(n — 2) -dimensional Liouville tori of the reduced system (7.6). 

So, in this approach, after solving the pendulum type system, the equations of Hcss- 
Appel'rot rigid body problem reduces to (dim SO (n — 2) + 1) differential equations of 
the reconstruction problem. 

Similarly as in subsection 5.1, it can be proved that the partial reduction of the rigid 
body system (4.10), (4.3), (4.15) is the pendulum system on the (n — l)-dimensional 
sphere S 1 ™ -1 (see [36]), given by the Hamiltonian function 

ho(x,p) = Jl + Jn (p,p) +Xn(x,7). 

Here the cotangent bundle of the sphere is realized as a submanifold (x, x) = 1, 
(x,p) = of R 2ll (x,p) and x, 7 € 5 n_1 represent the direction of the position of the mass 
center x an d the position of the vertical axes T in the space coordinates, respectively. 

7.3. Mishchenko Fomenko flows. Wc follow the notation of subsections 4.3 and 
5.4, where we take (M, G, K, P, N ) = (T*G, G a , G, g, T*(G/G a )). Here we consider the 
right G a -action and the left G-action on T*G. The momentum maps read $(<?,£) = 
pr 0a (£) and *&(g,£) — Ad s £, respectively. 

The metric K a ^, c is right G a -invariant and we can project it to the homogeneous 
space G/G a , that is to the adjoint orbit of a. Since we deal with the right action, the 
vertical distribution is left-invariant: V g = g ■ g a , while from the definition of K a ,b,c, the 
horizontal distribution is Ti g = g ■ d and the submersion metric does not depend on c. 
The submersion metric is exactly the metric K a ^ on 0{a) defined above. 

The Hamiltonian H a ^,c',D defines left-invariant metric on G that we shall denote 
by n a ,b,c',D- The Hamiltonian functions H a ^,c,D and H a ^,c coincides on the invariant 
manifold 

(7.7) (rG)„~,Gx5. 
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Therefore, the partial reduction of the geodesic flow of the metric K a ^,c,D is the 
geodesic flow of the metric K a ,b- The flow (7.3) is completely intcgrable for any e. 
Thus the geodesic flow of the perturbed Mishchcnko-Fomenko metric K a ,b,c,D satisfies 
Zhukovskii property. 

Remark 7.1. If a is a regular element of the Lie algebra g (G a is Abelian) then the 
invariant manifold (7.7) is almost everywhere foliated on invariant isotropic tori of the 
geodesic flow n a .b.c ,d and the motion over the tori is not quasi-periodic (see Example 
5.1). The same holds for the left G-reduced flow (4.30) restricted to (4.31). Namely, it 
can be proved that S C C and that C is complete at a generic point £ that belongs to 
(4.31) (e.g, see Theorem 3.2 in [12]). According Lemma 6.1, the perturbed flow (4.30) 
has the same foliation of (4.31) on invariant tori as the non-perturbed flow (4.25). Item 
(ii) of Theorem 6.1 then implies the integrability of geodesic flow of the metric K a ^. 

Remark 7.2. For a singular a G g, C + S restricted to (4.31) is not complete. The 
complete integrability of the geodesic flow of the normal metric and the magnetic geodesic 
flows (7.3) on the adjoint orbit 0{a) follows from the completeness of the commutative 
set Co induced from (4.27) within the algebra Cq(T*0(o)). The completeness is obtained 
by verifying the condition stated in Proposition 6.1 for the pencil of compatible Poisson 
structures within C£(T*0(a)), induced from (4.26) and (4.28), see [12, 58, 13]. Besides, 
item (hi) of Theorem 6.1 gives us an estimate of the dimension of invariant manifolds 
within (4.31) of the perturbed flow (4.30). 

Remark 7.3. The horizontal distributions H! and H (see subsection 5.5) of the metrics 
K a ,b,c,D and K a ,b,c are different 



be the integrals of the perturbed system (4.30) not only on the invariant manifold (4.31)) 
but on g as well. 

7.4. Singular Manakov flows. Similar perturbations as those of the Mishchcnko- 
Fomenko flows can be performed for other integrable Euler equations with symmetries. 
The natural candidate is the singular Manakov flow. Let 



be the orthogonal decomposition, where g = {X G so(n) | [X 7 a] = 0}. By M s and M a 
we denote the projections of M G so(n) with respect to (7.8). Further, let C : g — > g be 
an arbitrary positive definite operator. We take a and b such that the sectional operator 
A a ,b,C '■ so(n) — ► so(n) defined via 



is positive definite. Here ad and adt, are considered as invcrtible linear transformations 
from to [a, 0] C Sym(n). Let H a _b,c = A a ^,c(M)). We refer to Euler equations 





(7.9) 



{M* + M 3 ) = ad- 1 ad 6 (M 8 ) + C(M a ) 



(7.10) M — [M, O] , n - VH aAC {M) = A aAC (M), 
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as the singular Manakov Flow. The operator A a ,b.c satisfies Manakov condition [M, b] = 
[O, o], so we have the Lax representation with rational parameter A (sec Manakov [50]): 

(7.11) L(X) = [L(X),U(X)}, L(X) = M + Xa, U(X) = fl + Xb. 

In the case the eigenvalues of a are all distinct, i.e., = 0, Manakov proved that the 
solutions of the Euler equations (7.10) are expressible in terms of 9- functions by using the 
algebro-geometric integration procedure developed by Dubrovin in [24] (see [50]). The 
explicit verification that integrals arising from the Lax representation 

(7.12) C = {tr(M + Xa) k \ k = 1, 2, . . . , n, X e R}, 

form a complete commutative set on so(n) was given by Mishchenko and Fomenko [52]. 

Let us denote the set of linear functions on q by S. These are additional integrals 
in the case the eigenvalues of C are not all distinct and C is proportional to the identity 
operator, or more generally in the case C is an Adc-invariant, where G = SO(ki) x 
SO(k 2 ) x • • • x SO(k r ) C SO(n). The complete integrability of the system is proved 
by Bolsinov by using the pencil of compatible Poisson brackets given by the canonical 
Lie-Poisson bivector 

7Tl(6,6)|M = -(M, [6, 6]> 

and 

(see [10] and [70], pages 241-244). The another proof is given in [23]. Namely, we have 
{£,<S} so („) = and C + S is complete at a generic M e so(n). 

Now, the perturbation follows the perturbation of Mishchcnko-Fomenko flows: 

(7.13) M = [M, fi], Sl = VH aAC , D (M), 

where D : d — > g and H a ^ t c,D — H a ,b,C + (M s ,D(Mj,)) is positive definite (we do not 
suppose that C is Adc-invariant). The system (7.13) has the invariant manifold 

(7.14) d : M B = 0. 

Besides, the restriction of the system to (7.14) has the Manakov L-A pair (7.11) and 
integrals (7.12). 

Consider diagram (5.14), where we take 

(M,G,K,P,N ) = {T*SO(n),G, SO(n),so(n),T*(SO(n)/G)), 

with the right G-action and the left 5 l O(n)-action on T*SO(n). By the use of the map 
(6.4), the commutative set of function C induces a complete commutative set Co within 
Cf 0(n) {T*(SO{n)/G)) [23]. On the other side, if G is not commutative then C + S is 
not complete at 0. 

Proposition 7.1. If G is commutative, i.e., ai < 2, i = 1, . . . , r, then the set of function 
C + S is a complete set at a generic point Med. 

Proof. The proposition directly follows from item (iii) of Theorem 6.1. Alterna- 
tively, let L M = {V M tr(M + XA) k | k = 1, 2, . . . , n, A G R}. According to (9.3), C + S is 
complete at M if 

(7.15) (L M + &V 1 C L M + 9- 

The relation (7.15) is proved in [23] by using Theorem 9.2, namely by verifying that 
the dimension of the linear spaces (25) and (26) in [23] are equal to n. The dimension of 
the space (26) calculated in Lemma 2 [23] holds for a generic Me5. On the other side, 
the dimension of (25) is equal to n for elements in a generic position with the property 
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that M g is a regular element in d (sec [70], pages 234-237). Since g is commutative, the 
dimension of (25) will be n as required. □ 

Example 7.1. As an example, take n = 2r, a = (oti, a\, . . . , a r ,a r ). Then g = so(2) <g> 
so(2) ■■■(g) so(2) is the Cartan subalgebra. The set of integrals C + S is complete at J> 
so the invariant set d of the systems (7.10) and (7.13) is foliated on invariant tori. The 
matrix L(X) satisfies 

L12 = L2I = -^34 = ^43 = ■ • • = L 2r -l,2r = L 2r ,2r-1 = 0. 

In other words, the systems (7.10), restricted to (7.14), is an example of integrable iso- 
holomorhic system ([19], see Remark 8.2 given below). 

The bi-Hamiltonian formulation of singular Manakov flows can be performed by using 
the pencil II = {n\ lt \ 2 = Ai7Ti + A 2 7r 2 } given above (see [10, 70]). The singular brackets 
within the pencil II are proportional to ni t - ai , . . . , TTi,- ar and the linear function on 
so(ki) C g are among the Casimirs of the brackets wi- ai , i = 1, ■ ■ ■ ,r. So, if only one ki 
is greater of 1 (say k\ = k > 1, k 2 = • ■ • = k r = 1, r = n — k + 1, i.e., g = so(fc)), the 
perturbed singular Manakov flow (7.13) satisfies axiom (BP) with respect to the second 
Poisson structure 7Ti ; _ ai . The partially reduced system is completely integrable geodesic 
flow on the Stiefel variety SO(n)/SO(k). 

8. Integration of the magnetic pendulum on Gr + (4, 2) 

Let us consider closely the pendulum system given by the equations (7.6) in dimen- 
sion four with the magnetic term added. Since the orbit 0(x) — Gr+(4, 2) is defined 
with invariants, the cotangent bundle of the Grassmannian Gr+(4, 2) is given by the 
constraints: 

x 12 + x 13 + x 14 + x 23 + x 24 + ^34 = Xl2 + X34: 
^ £34X12 + X23X14 - 213X24 = X12X34, 

212P12 + X 13 pi 3 + X14P14 + X 23 P23 + 2J24P24 + X34P34 = 0, 
X34P12 + X23P14 ~ X 13 p 2 4 + P34^12 + ^23^14 - ^13^24 = 0. 

Introducing the magnetic momentum mapping (7.1), the equations 

x = (Ji + J n )[x, \p,x]\, 

p= (Ji + J n )\p, \p,x\] +e[x,p] - 7 + pr so(4)x 7 

become: 

m = [7, x] 

(8 2) 

x = (Ji + J n )[m,x]. 

The equations (8.2) are special case of the equations of the completely symmetric 
Lagrange bitop. The Lagrange bitop is defined in [18] and studied in details in [19]. 
Thus, the integration procedures given in [18, 19] can be applied to the considered 
system. We will present here both of them, the classical and the algebro-geometric 
integration procedures. 

After solving the system (8.2) one has m and x as known functions of time. In order 
to find p as a function of time, one needs to solve the equation 

m = [x,p] + ex 

in p. Let us denote 

too = m — ex = [x,p\. 
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Since (\p,x],so(4) x ) = ([x, so(4) x ],p) = we have too <= so (^)x- The operator ad^ : 
so(4) x i > 50(4)^ is bijective, thus a solution p = ad~ 1 (m ) is unique. 

8.1. Classical integration of the magnetic pendulum. Starting from a well- 
known decomposition so(4) = so(3) © so(3), let us introduce as in [19] 

TOl = i(TO + +TO_), TO 2 = i(TO+-TO_) 

where to + ,to_ are two three-dimensional vectors which correspond to four-dimensional 
matrix to^ according to 



(to+,to_) 



/ 





— m+ 


2 

TO+ 


— toL^ 


TOj_ 





— TO^ 


—mi 




2 


TO^ 





-TO 3 


V 


TO 1 ! 


2 

m_ 


m 3 _ 


/ 



(Similar decomposition can be performed for a;, 7). 
Equations (8.2) become 

(8.3) rhj = 2(7i x Xi), = 2(Ji + J„)(ra, x a;*), i = 1, 2, 

where 

71 = (0,0, --(712 +734)), 72 = (0,0, --(712 -734))- 

If we denote mi = (pi, to 2 = (p 2 , 92, r 2 ), then the first group of the equations 

(8.3) becomes 

Pi = -27(1)3^(1)2, P2 = -27(2)3^(2)2, 
Ql = 27(1)3^(1)1, 92 = 27(2)3^(2)1, 

fi - 0, r 2 = 0, 

where we denoted with x^j the j component of the vector Xi. 
The integrals of motion are for i — 1,2: 

Tt = III 

(J\ + J n ){p 2 i + qf) + 27^)32^)3 = hi 

PiX(i)i + qiX(i) 2 + nx (i)3 = f i3 = ex(i)3 

2 1 2 1 2 _ r 2 

x (i)l + x (i)2 + x (i)3 — Jii — X(i)3' 

The constants and fa are found from the conditions (8.1). 

Following [19] and introducing pi,(Ji, defined with pi — pi cos cr^, qi = pi sin cr^, we 

get 

Al + pfo 2 = 47( 2 4 )3(X( 4 )3 ~ 
pf&i = 27 (i ) 3 (ex(i)3 ~ fiiX(i)3)- 
It follows that Ui = satisfy equations 

= Pi(ui), 

where 

P;(m) = -( Ji + J n ) 2 u a + u 2 B t + uC t + A, « = 1, 2; 

and 

Bi = 2/ i2 (Ji + J n ) - fa(Ji + Jnf, 

°i = 4 7( 2 ) 3 X( 4 )3 - fi2 + (^1 + Jnf fil{filfi2 - 2eX( l )37(»)3), 

A = -(2e7 (i)3 X(i) 3 - /ii/i2) 2 , * = 1,2. 
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So, the integration of the system 

dui 



f dui =t f. 
J y/WM ' J 



du 2 



' J 

leads to the functions associated with the elliptic curves E\ , E 2 given with: 
(8-4) Ei =E i (J 1 , Jn,/»i,7(»)3)e.X(i)3)/»3) = y 2 = Pi(u). 

The equations (8.3) are very similar to those for the symmetric top, and they are spe- 
cial case of the equations of the Lagrange bitop (see [19]). From the equations (8.3) one 
concludes that the dynamics of the magnetic pendulum on Gr + (4, 2) splits on two inde- 
pendent systems on the sphere S 2 . This splitting corresponds to the fact that Gr + (4, 2) is 
a product of two spheres. Let us mention that a general Lagrange bitop is more complex 
since it doesn't split on two independent Lagrange tops. 

8.2. Algebro-geometric integration procedure of the magnetic pendulum. 

Algcbro-geometric integration completely follows paper [19] (see also [21]). 
The starting point in the integration is the following Lax representation: 

Proposition 8.1 ([19]). The equations (8.2) has the Lax representation: 

L(X) = [L(X),A(X)} 

where L(X) = X 2 c — Xm + x, A(X) = X-f — (Ji + J n )m and c = Ji+J 7 

Remark 8.1. The magnetic spherical pendulum on adjoint orbits (7.2) admits a similar 
Lax representation that provides a complete set of commuting integrals for a regular c E Q 
[14]. " 

We will change the coordinates in order to diagonalize the matrix j^rj 7- In this 
new basis the matrices L(X) have the form L(A) = U^ 1 L(X)U, 



L(X) 



where A i2 = A 2 ci 2 



( -«A 34 


03 - ifii 
\i0*3 + 01 

Xm.12 + £12, 





-ifo + 04 
01 + iPl 

A34 = A 2 c 3 4 
1 

x 3 = 7T 



-PI - %Pl 10 3 -P 4 \ 

-iPt-Pt -Pz + iPa 
-iAi2 
iAi2 



— A77134 + X34, and 
(X13 + ix 23 ) , 



(8.5) 



2/3 = 



2/4 



1 / ■ 

~2 ( m i3 + ^23) , 



1 



03 = %3 + Aj/3, 

04 = x 4 + Xy 4 , x 4 — ^ (xi4 + ix 24 ) , 

03 = x 3 + Ay 3 , 
01 = x 4 + Xy 4 , 

The spectral polynomial p(A, fi) = dct (^L(X) — /j, ■ lj has the form 

(8.6) p(X^) = ^ + P(X)^ 2 + [Q(X)] 2 , 
where 

(8.7) P(X) = A 2 2 



(mi4 + im 24 ) . 



A34 



We can rewrite it in terms of ; 



40 3 0* 3 + 40 4 0* 4 , Q(A) = A 12 A 34 + 21(0^04 - 0301)- 
and xa: 



P(X) = AX 4 - BX 3 + DX 2 -EX + F, Q(A) = GA 4 - HX 3 



IX 



JX + K. 
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Their coefficients 

^ = c 12 + c 34' 

B = 2c 34 m 34 + 2c 12 m 12 , 

D = mf 3 + mf 4 + to 23 + m\ 2 + m\ 4 + 2c 12 x 12 + 2c 34 x 34 , 

E = 2x 12 mi 2 + 2xi 3 TOi 3 + 2a;i 4 mi4 + 2a; 2 3^23 + 2:e 24 to 2 4 + 2X34TO34, 

F2 1 2 1 2 1 2 1 2 1 2 
= a; 12 -+- x 13 + x 14 + x 2 3 + a; 24 + a: 34 , 

G = C12C34, 

H = c 34 m, 12 + c 12 m 34 , 

I = c 34 x 12 + x 34 c 12 + TOi 2 m 34 + m 23 mi 4 - mi 3 m 24 
J = m 34 x 12 + m 12 x 34 + m 14 x 23 + m 23 x 14 - a;i 3 m 24 - x 24 m 13 , 
K = x 34 xi 2 + x 23 x 14 - x 13 x 24 

are integrals of the motion. From the constraints (8.1) one can calculate values of four 
integrals 

E = 2eF = 2e( X j 2 + X34), J = 2eK = 2e X i2X34. 
There is an involution a : (A,^) — > (A, — [i) on the curve T : p(A, ^1) — 0, which 
corresponds to the skew symmetry of the matrix L(\). Denote the factor-curve by Ti = 
T/a. 

Detailed analysis of algebro-geometric properties of the curves T, Ti one may find in 
[19]. 

We consider the next eigen-problem 

(J^+ A(X)^j tp k = 0, L(X)tp k = n k ip k , 

where ^fc are the eigenvectors with the eigenvalue fik- Then tpk(t,X) form 4x4 matrix 
with components tp k (t,\). Denote by corresponding inverse matrix. 
Let us introduce 

(there is no summation on k) or, in other words g(t) = V'fc(i) ®<p(t) k . Matrix g is of rank 
1, and we have 




We can consider vector-functions ipk(t, A) = {^\{t, A), ip k (t, A)) as one vector-function 
ip(t, (A, n)) — (tp 1 ^, (A, n)), ip 4 (t, (A, fi))) T on the curve T defined with ^(t, (A, fi k )) = 
ip k (t,X). The same we have for the matrix Lp\. The relations for the divisors of zeroes 
and poles of the functions ip l i ipi in the affine part of the curve T are: 

(8.9) (g}) a = d 3 (t)+d*(t)-D r -D> s , 

where D r is the ramification divisor over A plane(see [24] and [19]) and D s is divisor of 
singular points, D' s < D s . One can easily calculate deg_D r — 16, degZ) s = 8. 

The matrix elements g* (t, (A, )i k )) are meromorphic functions on the curve T. We 
need their asymptotics in the neighborhoods of the points P kl which cover the point 
A = 00. 

It was justified in [19] that from now on we may consider all the functions in this 
section as functions on the normalization T of the curve T. 
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Let us denote by dj and by d % the following divisors: 

d 1 =d 1 +P 2 , d 2 = d 2 +P u d 3 = d 3 + P i , d 4 = d 4 + P 3 , 
d l =d 1 +P 2l d 2 = d 2 + P 1 , d 3 =d 3 + P 4 , d A ^d i + P 3 . 
We have (see [19] for details): 
Proposition 8.2 ([19]). (i) The divisors of matrix elements of g are 
(g)) = <f + dj - D r + 2(P 1 +P 2 +P 3 + P^ -Pi- Pj 
(ii) The divisors di , d? are of the same degree 

deg di = deg d? = 5. 
Let us denote with $(i, A) the fundamental solution of 

*(i,A)=0, 



normalized with <j>(r) = 1. Then, if we introduce functions 

r(t, r, (a, /i fc )) = ^ $• (t, \)h s (T, (a, /i fc )) 

s 

where h s are the eigenvector of L(X) normalized by the condition J2 S (A, Mfc)) = 1) 
it follows that 

* (t, r, (A, w) ) - L A ) E^A) " E^A) ■ 

Proposition 8.3 ([19]). TTie functions satisfy the following properties 

(i) In the affine part of T the function ip l has 4 time dependent zeroes which belong 
to the divisor d l (t) defined by formula (8.9) , and 8 time independent poles, e.g. 

(^(t, t, (A, Mfe ))) a - d\t) - V, degf> = 8. 

(ii) At the points Pk, the functions ip l have essential singularities as follows: 

j>%T,(\,ii)) =exp [-(t-T)R k ]a%T,(\,n)) 
w/iere i?fc are given with 

-Ri = i y— (Ji + J n )m 34: j , R 2 = -R u R 3 = i (Ji + J n )m 12 j , R 4 = -R 3 

and a 1 are holomorphic in a neighborhood of Pk, 

&\t, t, (A, M )) = h\r, (A, fi)), &% r, P k ) = + vl(t)z + 0(z 2 ), 

with 

(8-10) v{ = 

C-ii Ckk 

We have 

Lemma 8.1 ([19]). The following relation takes place on the Jacobian Jac(T): 

A(d j (t) + ad j (t)) = A{d 3 (r) + ad? (r)) 
where A is the Abel map from the curve T to Jac(T). 



34 



DRAGOVIC, GAJIC, JOVANOVIC 



From the previous Lemma we see that the vectors A(d l (t)) belong to some translation 
of the Prym variety II = Prym(T | Ti). More details concerning the Prym varieties one 
can find in [57, 56, 29, 68, 69, 19]. 

It was shown in [19] that the Baker- Akhiezcr function '5 satisfies usual conditions 
of normalized (n=)4-point function on the curve of genus g = 5 with the divisor V of 
degree deg£> = g + n — 1 = 8, sec [26, 25]. By the general theory, it should determine 
all dynamics uniquely. 

Let us consider the differentials Q*. = gijdX, i,j — 1, ... ,4. 

It was proven by Dubrovin in the case of general position, that ttj is a meromorphic 
differential having poles at Pi and Pj, with residues v* and — v\ respectively. But here 
we have 

Proposition 8.4 ([19]). The four differentials Q\, f2|, f2| are holomorphic during 
the whole evolution. 

Remark 8.2. The proof was based on the fact that 
(8.11) v\ = v \=^ = vi = Q, 

which is consequence of condition L12 = L21 = £34 = L43 = 0. It was the reason that 
the notion of isoholomorphic systems has been introduced in [19] to describe such class 
of integrable systems. 

Let us recall ([19]) the general formulae for v. 

(812) XAA{P,)-A{P l)+ tU + Za ) 

where U = Y,x {k) U^ is certain linear combination of b periods of the differentials 
of the second kind ftp), which have pole of order two at Pi; A, are nonzero scalars, and 



[v\{A{Pi ~ Pj)) 



(-d u(i) 6[v}(o))V*(-d uW e[v}m i/2 y 

(Here v is an arbitrary odd non-degenerate characteristics.) Thus, from (8.12), it follows 

Proposition 8.5 ([19]). Holomorphicity of some of the differentials flj implies that the 
theta divisor of the spectral curve contains some tori. 

In a case of spectral curve which is a double unramified covering 

7r:f-r i; 

with .g(ri) = g, g(T) — 2g — 1, as it is satisfied for the Lagrange bitop, it is really 
satisfied that the theta divisor contains a torus, see [57]. Following [57] and [19], let us 
denote by LT~ the set 

r = {ie Pic 2g - 2 T I NmL = K ru h°(L) is odd} , 

where Kr 1 is the canonical class of the curve Ti and Nm : PicT — ► PzaTi is the norm 
map, sec [57, 69] for details. For us, it is crucial that II~ is a translate of the Prym 
variety II and that Mumford's relation ([57], p. 241-242) holds 

(8.13) n- C e f . 
Let us denote 

(8.14) U = i( X 3aU^ - X34 [/(2) + Xl2 [/(3) _ Xl2 [/W), 
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where f/W i s the vector of b periods of the differential of the second kind Op , which is 
normalized by the condition that a periods are zero. We suppose here that the cycles 
a, b on the curve T and a, b on Ti are chosen to correspond to the involution a and the 
projection ir, see [69]: 

7r(a ) = a ; 7r(6 ) = 2fo , v(a k ) = a k+2 , k = 1,2. 

The basis of normalized holomorphic differentials [uo, . . . , W5] on f and [v , v\, V2] on Ti 
are chosen such that 

TT*(v a )=U , 7T*(Ui) = Vi + a(Vi) — Vi + v i+2 , i = 1,2. 

Now we have 

Theorem 8.1 ([19]). (i) If the vector zq in (8.12) corresponds to the translation of the 
Prym variety II to IT 7 and the vector U is defined by (8.14) than the conditions (8.11) 
are satisfied. 

(ii) The explicit formula for z n is 

(8.15) z = ^(roo,roi,fo2,Toi,f 2), m = [ Ui, i = 0,1,2. 

The evolution on the Jacobian of the spectral curve, as we considered Jac(T) gives the 
possibility to reconstruct the evolution of the Lax matrix L(X) only up to the conjugation 
by diagonal matrices. As it was explained in [19], for complete integration one has to pass 
to the generalized Jacobian, obtained by gluing together the infinite points. Those points 
are Pi, P2, P3, P4 and corresponding Jacobian will be denoted as Jac(T\ {Pi, P2, P3, Pi})- 

It can be understood as a set of classes of relative equivalence among the divisors 
on r of certain degree. Two divisors of the same degree D\ and D2 are called equivalent 
relative to the points P\, P2, P3, P4, if there exists a function / meromorphic on T such 
that (/) = Di - D 2 and /(Pi) = /(P 2 ) = /(P 3 ) = /(P 4 ). 

The generalized Abel map is defined with 

A(P) = (A(P),Ai(P),...,A 4 (P)), A,(P)-cxp/ Op iQo ,i = l,...,4, 

J Pa 

and A(P) is the standard Abel map. Here fip,^ denotes the normalized differential of 
the third kind, with poles at Pi and at arbitrary fixed point Qo- 

We will use the generalized Abel theorem as it was formulated in [19]. The gener- 
alized Jacobi inverse problem can be formulated as the question of finding, for given z, 
points Qi, . . . , Qs such that 

8 4 

1 2 

8 rQs 

Xj = cexp V / ttpjQo + Kj,j = 1, ...4, 
s =i J p« 

where K is the Riemann constant and the constants Kj depend on the curve T, the points 
Pi, P 2 , P3, P4 and the choice of local parameters around them. 

We will denote by Q s the points which belong to the divisor T> from the Proposition 
8.3, and by E the prime form from [29]. Then we have 
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Proposition 8.6 ([19]). The scalars Xj from the formula (8.12) are given with 

8 rQs 

Xj = X°j cxp E ix ik) 7* , A° = c cxp E / to PiQo + Kj , 

k^j 8=1 

where the vector x — {x^ x \...,x^) denotes ^734, —734, 712, —712) <wid 

(kj 1 is a local parameter around Pj.) 

To give the formulae for the Baker- Akhiezer function, we need some notations. Let 

where 

7^= / m^i, 

and 7™ is defined by the expansion 
p 



Denote 



/ n<£ =-k m + ^ + o(k- 1 ), p^p m . 

JPo 

Finally, one can state 

Proposition 8.7 ([19]). The Baker- Akhiezer function is given by 

x o 8(A(P-P 3 )-z ) 

r(x,P) = efS(x,P) — 4 J JImpU)-^ 3 = 1 >->*> 

2^k=l A k e{P,P k ) 

where zq is given by (8.15). 

9. Appendix: Basic notions of the Hamiltonian systems 

9.1. Hamiltonian systems. Let (P, {■,■}) De a Poisson manifold and it be the 
associated bivector field on P 

{fi,f2}(x)=n x (df 1 (x),df 2 (x)) = J2^ i ^^ 1 



dxi dxj 

1,3 J 



If 7r is non-degenerate, then the two-form w = ^LOijdxi A dxj (uijn 3 ^ — 5f) is a 
symplectic form and (P, u>) is called a symplcctic manifold. 
The equations: 

(9.1) x = X h (x) <i=> f = {f,h}, /e<7°°(P) 

are called Hamiltonian equations with the Hamiltonian function /i and A£ = n l3 dh/dxj 
is the corresponding Hamiltonian vector field. 

A function / is an integral of the system (constant along trajectories of (9.1)) if and 
only if it commutes with h: {h, /} = 0. From the Jacobi identity the Poisson bracket of 
two integrals is again the integral, so we can consider a Poisson subalgcbra T C C°° (P) 
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of integrals (or a collection of integrals closed under the Poisson bracket). Consider the 
linear spaces 

(9.2) F x = {df(x) |/€f}c T* X P 

and suppose that we can find / functionally independent functions fi, ■ ■ ■ , fi G T whose 
differentials span F x almost everywhere on M and that the corank of the matrix {/j, fj} 
is equal to some constant fc, i.e., dim ker tt x \ p x = k. The numbers / and k are called 
differential dimension and differential index of T and they are denoted by ddim T and 
dind J 7 , respectively. 

We say that T is complete at x if the space F x given by (9.2) is coisotropic: 

(9.3) f: c F x . 

Here F x is skew-orthogonal complement of F x with respect to 7r: 

F: = {Z£T*P\tt x (F x ,Z) = 0}. 

The set T is complete if it is complete at a generic point x G P. In this case 
F x — ker it x \f x and dind T = dimF x , for a generic x G P. Equivalently, T is called 
complete if (see [10, 11, 74]): 

ddim T + dind T = dimP + corank {•, •}. 

The Hamiltonian system on (9.1) is completely integrable (in noncommutative sense) 
if it possesses a complete set of first integrals T. Then (under compactness condition) P 
is almost everywhere foliated by (dind T — corank {•, -})-dimensional invariant tori. As 
in the Liouville- Arnol'd theorem [6] , the Hamiltonian flow restricted to regular invariant 
tori is quasi-periodic (see Nekhoroshev [59], Mishchenko and Fomcnko [53] and Zung 
[74]). 

9.2. Natural mechanical systems. The basic examples of Hamiltonian systems 
are natural mechanical systems (Q,n,v), where Q is a configuration space, n is a Rie- 
mannian metric on Q and v : Q — > M is a potential function. Let q = (q 1 ,. . . ,q n ) be 
local coordinates on Q. The motion of the system is described by the Euler-Lagrangc 
equations 

,„ d dl dl 

(9.4) — — = — , i = l, ...,n, 
v ' dtdq 1 dq*' 

where the Lagrangian is l(q, q) = \{n q q, q) - v(q) = \ Kijtfq 3 - v(q). 

Equivalently, we can pass from velocities cf to the momenta pj by using the standard 
Legendre transformation pj — n^q 1 . Then in the coordinates q l ,Pi of the cotangent 
bundle T*Q the equations of motion read: 

dq l dh dpi dh . 

(9 - 5) -M=dp~^ ^ = ~W 

where the Hamiltonian h is the sum of the kinetic and potential energy of the system 
h(q,p) = \ J2ij K ^PiPj + v (q)- Here K y are the coefficients of the tensor inverse to the 
metric. 

This system of equations is Hamiltonian on T*Q endowed with the canonical sym- 
plectic form w = Ym=i d-Pi A dq % . The corresponding canonical Poisson bracket is given 
by 



(9-6) {/,<?} = £ 



df dg _ dg df 
dq 1 dpi dq l dpi 
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Let e be a real parameter (a "coupling" constant). The motion of the particle 
under the influence of the additional magnetic field given by a closed 2-form eft = 
J2i<i<j< n e Fij{l)dq l A dqi , is described by the following equations: 

dcf__ dh_ dp 1 _ _dh_ ^ dH 
{ ' dt ~ dpS dt dqi +6 f^ l3 dpi 

The equations (9.7) are Hamiltonian with respect to the symplectic form u) + ep*fi, 
where p : T*Q — > Q is the natural projection. Namely, the new Poisson bracket is given 

by 

(9.8) {/.^ = {/ lS } + e^,|| 

and the Hamiltonian equations / = {/, h} e read (9.7). 

9.3. Hamiltonian G-actions. Let a connected Lie group G act on 2n-dimensional 
connected symplectic manifold (M, co). The action is Hamiltonian if G acts on M by 
symplectomorphisms and there is a well-defined momentum mapping: 

(9.9) $ : M -> fl * 

(g* is a dual space of the Lie algebra g) such that one-parameter subgroups of symplecto- 
morphisms are generated by the Hamiltonian vector fields of functions 4>^(y) = (3>(y),£), 
£ G g and (p^,^] = {4>Zn 4>& }■ Then $ is equivariant with respect to the given action of 
G on M and the co-adjoint action of G on g*: • x) = Ad*($(x)). In particular, if rj 
belongs to $(M), then the co-adjoint orbit 0{rj) belongs to $(M) as well. 

The mapping / i— > / o $ is a morphism of Poisson structures: {/i o / 2 o 3>}(a;) = 
{/l, /2}g* (n), = ${x), where {•, -} B * is the Lie-Poisson bracket on g*: 

{/i,/2} fl .(»?) - (r?,[#i(r?),rf/ 2 (r?)]), /i,/ 2 : 0* - M- 

Thus, $*C°°(g*) is closed under the Poisson bracket. Since G acts in a Hamiltonian 
way, the set of G-invariant functions Cq(M) in G°°(M) is closed under the Poisson 
bracket as well. Also {<S>*C°°(t*),C%>(M)} = (the Nocthcr theorem). 

Suppose the group G is compact. Then we have 

Theorem 9.1 ([11]). (i) The algebra of functions $*C°°(g*) + Qf(M) is complete: 

ddim ($*G°°(g*) + Ca (M)) + dind (**G°°(fl*) + C£?(M)) = dimM. 

(ii) Suppose A C G°°(fl*) is a involutive set of functions, complete on a generic 
coadjoint orbit 0{rf) C <&(Af) and B is a complete commutative subset of Cq(M): 

ddim B = dind B = ^ (ddim Cg'(Af) + dind Cg'(M)) . 

T/ien $M + C£?(M) and $*G°°(g*) + B are complete sets on M, while $*A + B is a 
complete commutative set on M. 

9.4. Symplectic reductions. Let G be a Lie group with a free and proper Hamil- 
tonian action on a symplectic manifold (M, ui) with the momentum mapping (9.9). As- 
sume that r\ is a regular value of $, so that M v — ^ 1 (n) and Mo, = ^~ 1 (O v ) are 
smooth manifolds. Here O n = G/G v is the coadjoint orbit of n. The manifolds M n 
and Mo v are G^-invariant and G-invariant, respectively. There is a unique symplectic 
structure u v on N v = M^jGr, = ^ 1 {0 71 )/G satisfying 



SYSTEMS OF HESS-APPEL'ROT TYPE AND ZHUKOVSKII PROPERTY 



39 



where n ri : M v — > N v is the natural projection (Marsden and Weinstein [47]). According 
to Noether's theorem, if h is a G-invariant function, then the momentum mapping $ is an 
integral of the Hamiltonian system x = Xh(x). In addition, the restriction of Xh to the 
invariant submanifold M v projects to the Hamiltonian vector field X^ v on the reduced 
space N v with h n defined by /i|m„ = = h n o ir v . 

9.5. Cotangent bundle reductions. As the important example, consider the tr- 
action on the configuration space Q. The action can be naturally extended to the Hamil- 
tonian action on (T*Q,u>): g • (q,p) = (g ■ q, with the momentum mapping $ 
given by 

(9.10) ($( g ,p), = (p, 

where t; q is the vector given by the action of one-parameter subgroup exp(t£) at q [42]. 

Now, let G be a connected Lie group acting freely and properly on Q and 7r : Q — > 
Q/G be the canonical projection. Then is the regular value of the cotangent bundle 
momentum mapping (9.10) and the reduced space ($ _1 (0)/G, uj ) is symplectomorphic 
to T*(Q/G). 

Suppose (Q, k, v) is a G-invariant natural mechanical system. That is G acts by 
isometries and the potential is the pull back of the potential V defined on Q/G. The 
metric k induce the submersion metric on Q/G (e.g., see [7]). The reduced system, for 
a zero value of the momentum mapping, is the natural mechanical system (Q/G, K,V). 
For Abelian groups this is the classical method of Routh for eliminating cyclic coordinates 
[67]. Within Lagrangian formalism the non- Abelian construction for the zero level-set of 
and for the other values of the momentum mapping are given in [5] and [49] , respectively. 

9.6. Compatible Poisson brackets. Let {•, -}i and {•, -} 2 be compatible Poisson 
structures on a manifold P. In other words, each linear combination Ai{-, -}i + A 2 {-, -} 2 
with constant coefficients is again a Poisson structure (e.g., see [10, 70, 72, 63] and 
references there in). Let tti and 7r 2 be the associated bivector fields and let 

n = {7TAi,A 2 | Ai, A 2 e M, A 2 + X% ^ 0}, 7r Al ,A 2 = Aitti + A 2 7r 2 . 

In what follows we shall suppose that all functions are defined on some open set U, 
x G U. By r denote the corank of a generic bracket (or bivector) in n at x. For each 
bracket in n of corank r, we consider the set of its Casimir functions at x. Let C be the 
union of these sets. Then C is involutive set with respect to every Poisson bracket from 
n. Let C x denote the linear subspace of T*P generated by the differentials of functions 
from C. It is clear that C x is spanned by the kernels ker7r(x), it 6 n, corank ir(x) = r. 

Together with II, consider its natural complexification II C = {7Ta 1 ,a 2 = 7TiAi +7r 2 A 2 , 
Ai,A 2 e C, |Ai| 2 + |A 2 | 2 7^ 0}. Here, we consider ir\ li x 2 as a complex valued skew- 
symmetric bilinear form on the complexification of the co-tangent space (T*P) C . There 
are only finite number of the non-proportional singular structures 7r A i i A i , . . . , tt x p ^ G n c 
with a corank greater then r at x. With the above notation, we can state the following 
remarkable result: 

Theorem 9.2 (Bolsinov [10]). (i) G£ does not depend on the choice it G n. 

(ii) (G£) c DGC + ker7r MiA i(a ; ) + --- + ker7r A P ;A P (x) . 

(iii) (G£) c = G^r + ker7r A i )A i(or) H h ker7r A P !A p (;r) if and only if 

(9.11) dim c K Ai;A , - r, ^ A « , A » - kev ^ C (7;*G) C , i = l,...,p, 
where ttq G n is of the maximal rank at x. 
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As a corollary, an important completeness condition is formulated in [10]: 

Theorem 9.3 ([10]). Let n £ II and corank ir(x) = r. Then C is a complete commutative 
set at x G P if and only if corank 7r'(x) = r for all tt' <E n c ; n' ^ \tt, A £ C. 

By the use of Theorem 9.2 one can also formulate conditions for non-commutative 
integrability in the case that some of the brackets in II are not of the maximal rank [10] 
(see also [63]). 
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